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We demonstrate that the skeleton of the Fermi surface S^-a pertaining to a uniform metallic 
ground state (corresponding to fermions with spin index a) is determined by the Hartree-Fock 
contribution S"^(fc) to the dynamic self-energy T.a(k;e). That is to say, in order for k G S^-a, it is 
necessary (but for anisotropic ground states in general not sufficient) that the following equation be 
satisfied: 



where stands for the underlying non-interacting energy dispersion and Ef for the exact interacting 
Fermi energy. The Fermi surface <Sf;ct consists of the set of k points which in addition to satisfying 
the above equation fulfil 



The set of k points which satisfy the first of the above two equations but fail to satisfy the second 
constitute the pseudogap region of the putative Fermi surface of the interacting system. We consider 
the behaviour of the ground-state momentum-distribution function ndk) for k in the vicinity of 
<Sf;o- and show that, whereas for the uniform metallic ground states of the conventional single-band 
Hubbard Hamiltonian, described in terms of an on-site interaction, ncT{k^.a-) > | and r\a{k^.„) < i 
(here k^.^ and kf.^ denote vectors infinitesimally close to <Sf;ct, located respectively inside and outside 
the underlying Fermi sea), for interactions of non-zero range these inequalities can be violated 
(without thereby contravening the stability condition ncr{k^.„) > na{k^.,j)). This aspect is borne 
out by the r\a{k) pertaining to the normal states of for instance liquid "^Ife (corresponding to a 
range of applied pressure) as determined by means of quantum Monte Carlo calculations. We 
further demonstrate that for Fermi- liquid metallic states of fermions interacting through interaction 
potentials of non-zero range (e.g. the Coulomb potential), the zero-temperature limit of ncr{k) does 
not need to be equal to i for k G Sf-o-', this in strict contrast with the uniform Fermi-liquid metallic 
states of the single-band Hubbard Hamiltonian (if such states at all exist). This aspect should be 
taken into account while analyzing the no-(fe) deduced from the angle- resolved photoemission data 
concerning real materials. We discuss, in the light of the findings of the present work, the growing 
experimental evidence with regard to the 'frustration' of the kinetic energy of the charge carriers in 
the normal states of the copper-oxide-based high-temperature superconducting compounds. 



Ek + hYfJ{k) = e, 



S'a(fc) =0, 



where 
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§ 1. INTRODUCTION 

Theoretical investigation of anisotropic interacting 
metallic systems is specifically demanding through the 
possibility, even for isotropic interaction potentials, of the 
deformation of the Fermi surface in these systems with 
respect to the underlying non-interacting Fermi surface 
[1-3] . For recent approaches to this problem we refer the 
reader to [4, .5] and the references herein. 

Building on recent work [6] concerning the uniform 
metallic ground states (GSs) of the conventional single- 
band Hubbard Hamiltonian [7 10] for arbitrary spatial 
dimensions d, in this work we consider similar states 
for arbitrary two-body interaction potentials (provided 
only that these possess Fourier transform) correspond- 
ing to systems of fermions whose non-interacting energy 
dispersion Sk can be both boimded (as is the case for 
the single-band Hubbard Hamiltonian) and unbounded 
from above (as in the case of the free-fermion model 
where Sk = ^^||fc|P/[2m], in which m stands for the 
fermion bare mass). Our results, of which a number we 
have presented in the abstract, make explicit some of 
the fundamental differences between systems in which 
the interaction potential v{r — r') is of contact type, 
of strictly bounded and of unbounded range. Follow- 
ing the general considerations leading to the equations 



defining sSpicr, presented above, in this paper we deduce 
the most general expression for ^^{k) in the close vicin- 
ity of <Sf;ct and for specific cases deal with this expres- 
sion in considerable detail. Amongst others, we show 
that, in contrast with the case of the conventional single- 
band Hubbard Hamiltonian, n<j(fe) can become smaller 
(greater) than ^ for fc inside (outside) the Fermi sea 
and inflnitesimally close to S^-^ai without infringing on 
the condition n^{k~.^) — n^{k'^.^) > implied by the as- 
sumption of the stability of the GS of the system under 
consideration. This aspect is significant in that in cases 
where ncrikp.^) < 5, by continuity n„{k) is expected to 
be more suppressed for fc inside the Fermi sea (at least 
in the vicinity of Sj,.^„) in comparison with cases where 
n<^(fe~^) > |, implying thus an increase in the contri- 
bution of the correlated kinetic energy g.£fcncr(fc) to 
the GS total energy of the system under consideration 
(see later). In this connection, note that ^ka'^'^i^) 
equal to the total number A'^ of particles in the GS and 
therefore independent of the magnitude of the coupling 
constant of interaction. 

For a number of reasons of contemporary interest, a 
reliable description of na{k) seems timely and desirable. 
Foremost amongst these are investigations regarding the 
determination of the mechanism of superconductivity in 
the copper-oxide-based high-temperature superconduct- 
ing materials (hereafter referred to as cuprates or cuprate 
compounds); for recent reviews see [11 13]. Both the ear- 
lier angle-resolved photoemission results [14] (see in par- 
ticular Fig. 3(F) herein) and the recent in-plane optical- 
conductivity data concerning Bi2Sr2CaCu208+5 [15,16] 
(see also [17]) reveal decrease in the in- plane kinetic 
energy of the system upon entering into the supercon- 
ducting state, by an amount commensurate with the ex- 
pected superconducting condensation energy. ^ A simi- 
lar trend has been observed in the measured out-of-plane 
(i.e. the c-axis) optical conductivities of Tl2Ba2Cu06+a;, 
La2-a:Sra;Cu04, and YBa2Cu306.6 [18]. ^ Assuming the 
sufficiency of the BCS variational wavefunction to de- 



^ This aspect has been explicitly quantified in [15-17]. 

^ The 'interlayer tunnelling' theory by Anderson and co- 
workers [19,20] (see also [21]) ascribes the superconducting 
condensation energy wholly to the reduction in the kinetic en- 
ergy of the cuprate compounds in consequence of the possibil- 
ity of coherent intcr-laycr tunnelling of the paired electrons in 
the superconducting states, a process deemed infcasible in the 
normal states of these materials (owing to the incoherent na- 
ture of the underlying single-particle excitations). Although, 
among other things (see [22] and [20]), the measured c-axis 
optical conductivity in Tl2Ba2Cu06+a;, La2-iSra;Cu04, and 
YBa2Cu306.6 [18] may be interpreted within the framework 
of the inter-layer tunnelling theory [18] (see however [23,24]), 
some experiments [25,26] appear [27] (see also [28]) not to 
support this theory; see however [29], [30] (and [31]). The 
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scribe the superconducting states of the cuprate com- 
pounds (an assumption that has been widely antici- 
pated), ^ one readily observes that n^{k^.^) < ^, with 
fepio- G >5p;o., in combination with the assumption of con- 
tinuity of na{k) for k inside the Fermi sea and in the 
vicinity of S^-cr, is tantamount to an excess kinetic energy 
in the normal GS of the imdcrlying system in compari- 
son with the kinetic energy in the paired state (see the 
previous paragraph). 

The above-indicated experimental results on the one 
hand and the feasibility of ncr{k~.^) < ^ for the metal- 
lic states of fermions interacting through potentials of at 
least finite range on the other clearly signify the likeli- 
hood of the existence of a crucial relationship between 
the latter potentials and the unconventional properties 
of the metallic states of the cuprate compoimds (see also 
[37,38]). In this connection it is relevant to mention 
that it has been rigorously shown [39] that for = 2 
the conventional single-band Hubbard Hamiltonian does 
not support dx2-y2 long-range pairing order at non-zero 
temperatures, ^ an aspect strongly anticipated by earlier 
Monte Carlo results (which incidentally also do not sup- 
port extended s-wave pairing) on finite systems [44-48]. 
Further, the view is increasingly gaining prominence that 
the long range of the Coulomb potential, such as taken 
account of by means of a variety of 'extended' Hubbard 
(and associated t-J) Hamiltonians [49-65], plays a cru- 
cial role in stabilizing one from amongst a number of 
competing orders in the cuprate compounds over some 



significance to pairing of the in-plane kinetic-energy frustra- 
tion in the normal states of the cuprate compounds has been 
put forward and established by the model considerations of 
Hirsch and Marsiglio [32,60] and Hirsch [33,52] (sec also [34]). 

^ According to this assumption, the zero-temperature BCS 
occupation fraction = 1 — u^, (see [35], Chapter 3 and in 
particular Fig. 3.1) plays the same role as nCT(fe) does in rela- 
tion to the underlying normal GS. For completeness, within 
the framework of the conventional BCS theory, and at zero 
temperature, the na{k) pertaining to the normal GS coin- 
cides with the unit-step function, equal to unity for k inside 
the Fermi sea and to zero for fc outside. Consequently, at 
zero temperature the BCS state accommodates a higher ki- 
netic energy than the underlying uncorrelated normal GS (for 
a quantitative consideration see § 3.4 in [35]). Note in pass- 
ing that Chester's [36] findings do not apply to a pairing, or 
reduced, Hamiltonian such as the one employed in the con- 
ventional BCS theory. 

* The superconducting order parameter in the cuprate com- 
pounds is generally considered to have the A^2_y2 symme- 
try [11-13]. Note in passing that in practice the inter-planar 
coupling of the electronic degrees of freedom is capable of 
stabilizing a d-wave superconducting state at finite tempera- 
tures, such as obtained within the framework of the so-called 
'cluster dynamical mean-field' approximation [40,41] (see also 
[42,43]). 



relevant range of parameters. Of these, one concerns the 
periodic (but incommensurate with respect to the under- 
lying lattice) ordering of spin and charge (also referred to 
as spin-charge stripes; for general reviews sec [66,67] and 
the references herein), which appears to compete against 
macroscopic phase separation [68]. This statement re- 
quires some qualification. Macroscopic phase separa- 
tion is believed not to concern the Hubbard Hamiltonian 
[69,70] , but solely the t-J Hamiltonian. However, even for 
the t-J Hamiltonian the actual situation is not firmly es- 
tablished; on the one hand, Green-function Monte Carlo 
results (concerning finite systems) have been interpreted 
as implying phase separation at hole dopings and inter- 
action strengths appropriate to the cuprate superconduc- 
tors [71] , while, on the other hand, subsequent numerical 
calculations [70] based on 'density-matrix rcnormaliza- 
tion group' (DMRG) method performed on larger sys- 
tems have upheld the viewpoint that the 'stripes' state 
would have a lower energy than the phase-separated 
state. In spite of these opposing observations, it is gen- 
erally believed that the sensitivity of the results corre- 
sponding to the t-J Hamiltonian on such details as, for 
instance, the shape and size of the clusters on which 
computations are carried out does not warrant an a pri- 
ori neglect of the long range of the Coulomb interaction 
[61,63,64] (see also [72]). Later in this paper we shall 
discuss and document the significance of a non-contact- 
type interaction in accurately describing the normal liq- 
uid state of ^He. 

Using the experimental angle-resolved photoemission 
data it is possible to determine n(j(fc) for k in the vicinity 
of iSf:<j, as has been reported in [73] for Bi2Sr2CaCu208-i-5 
(Bi2212) and BiaSraCuiOe+a (Bi2201), in both the over- 
doped and the underdoped regime. In [73] the Fermi sur- 
face has been asserted as being the locus of the k points 
at which, firstly, the (experimental) single-particle spec- 
tral function A^{k; e) is peaked for e = and. secondly, 
r\a-{k) takes the value i (see [74]); ' in [73] the abso- 
lute value of n„(k) is determined by equating with unity 
the maximum value achieved by the experimental n^rik), 
Vfc, a procedure which evidently adds uncertainty to the 
experimentally determined r\„{k). As demonstrated in 
[6], and in view of the findings in the present work, the 



^ Explicitly, in [73] the authors introduced three criteria to 
be obeyed "at a true FS [Fermi surface] crossing", of which 
the second reads: "(2) n{k) [r\a{k)] should be near 50% of its 
maximal value, or equivalently the maximal gradient point" . 
Elsewhere these workers stated: "a FS [Fermi surface] crossing 
should occur when n{k) loses about half of its maximum value 
(excluding the background)" . As we indicate in the main text, 
these workers fixed the scale of the experimentally determined 
nCT(fe) by identifying with unity the maximum value of the 
latter function over all the k points probed experimentally. 
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requirement r\a{k) = ^ for fc £ 5f;ct is not justified in gen- 
eral; ^ firstly, within the framework of the conventional 
single-band Hubbard model, the latter is solely a neces- 
sary (and not sufficient) condition for metallic states to 
be Fermi-liquid (FL) and, secondly, for metallic states of 
fermions interacting through potentials of non-zero range 
the condition n^ik) = i for fc G S^-a is not demanded 
even for FL metallic states. In the light of these and 
neglecting for the moment the above-mentioned uncer- 
tainty concerning the absolute values of the experimen- 
tally determined ri^{k), the results depicted in Fig. 2 of 
[73] have far-reaching implications. For instance, these 
results clearly show that a contact-type interaction can- 
not be capable of reliably describing Bi2212 and Bi2201 
in the over-doped and the underdoped regimes (at no 
peak position of Aa{k;eY) does r\a-{k) acquire a value 
reasonably close to a number of these results further 
show that the pertinent single-particle spectral functions 
Acr{k;e), corresponding to the systems considered, arc 
most probably free from quasi-particle ^-function contri- 
butions (establishing the breakdown of the FL state) at 
regions of Sp-a neighbouring the M points of the underly- 
ing Brillouin zones (compare for instance Figs. 2(g) and 
2(i) in [73] with the results in Eq. (109) below). 

The organization of this paper is as follows. In § 2 
we introduce the Hamiltonian on which the considera- 
tions in this paper are based; this Hamiltonian, which 
is expressed in terms of a single spin-degencratc energy 
dispersion £fc, is defined on the continuum; however, 
through some minor adjustments, which we explicitly 
specify, this reduces to the conventional single-band Hub- 
bard Hamiltonian, defined on a lattice. In § 3 we present 
the details underlying the equations from which S^-.a per- 
taining to the uniform GSs of the Hamiltonian introduced 
in § 2 is determined. In § 4 wc expose an intimate rela- 
tionship between S^-a and T,^^{k) mediated through the 
single-particle density matrix Qai^: f') which for uniform 
GSs, where the latter function depends on r — r', coin- 
cides with the Fourier transform of the GS momentum 
distribution function n^ik): the information with regard 
to <Sf;o-, through the singular nature of nij{k) for k G S-p-a^ 
is shown to be already fully reflected in the leading term 
in the asymptotic scries of £Io-(t', »"') for ||r — r'|| oo. In 
§ 5 we examine the equation for iSF:cr in terms of SJJ^(fc) 
by applying this equation to an isotropic metallic GS 
for which iSp;^ is rotationally invariant and is fully deter- 
mined by the concentration of the fermions in the system. 
By doing so we deduce expressions for the exchange po- 



® For cases in which nc{k) is discontinuous at fe € iSf;ct, in 
the present context the value assigned to no-(fc) for k € Sy;<t 
is the arithmetic moan of n^ik') for k' approaching fe from 
inside and outside the Fermi sea. This mean value coincides 
with the zero- temperature limit of n^ik) for fc € iSf;^ [6]. 



tential /ix:cr and correlation potential ^c-.a, as encountered 
within the framework of the GS density-functional theory 
[75-77]. Our expression for /Xx;<7 identically reproduces 
the exact /ix:cr corresponding to the Coulomb-interacting 
uniform-electron-gas system. Making use of n^^*(fc), the 
random-phase approximation (RPA) to ncr(fe), which for 
systems of Coiilomb-intcracting fermions asymptotically 
approaches the exact n(,{k) in the high-density regime, 
we calculate ^j,c;a which in the latter regime up to an 
unimportant additive constant reproduces the ^c^a de- 
duced from quantum Monte-Carlo calculations. In § 6 we 
consider the behaviour of ^^{k) for fe in the vicinity of 
5F:cr, both inside and outside the Fermi sea. Rather than 
being exhaustive, in this Section we attempt to contrast 
the important consequences that the range of interaction 
can have on the behaviour of n(j(fc) in the mentioned re- 
gion of the fc space. In § 7 we present a summary of 
the main aspects of this paper. Here we also briefly re- 
late our theoretical findings with some available results 
concerning correlated fermion systems, specifically liquid 
^He in the normal state. Finally, in appendix A which 
follows the main body of this paper we present an anal- 
ysis which exposes the nature of inaccuracy in (fe) as 
calculated on the basis of a flnite-order perturbation se- 
ries for the self-energy operator. For reasons that will be 
clarified in appendix A, the defining expression for S^-^a 
as introduced in this paper (i.e. Eq. (53) below) more 
reliably reproduces iSf;^ than the conventional defining 
expression (i.e. Eq. (45)). 

§ 2. PRELIMINARIES 

In a recent work [6] wc have in some detail consid- 
ered the uniform metallic GSs of the single-band Hub- 
bard Hamiltonian. Here we extend our investigations by 
considering uniform metallic GSs of the following Hamil- 
tonian 

ii = i?o + gHi, (1) 

where 

Ho = ^ekcl.„Ck-a, (2) 

^1 = ^ 13 13 4+9;aCp-q;a'Cp;a'Cfe;a- (3) 

<T,cr' k,p,q 

In Eq. (1), g (which has the dimension of energy) stands 
for the coupling-constant of interaction, Sk for a spin- 
degenerate single-particle energy dispersion (which may 
be equated with ?i^||fe|p/[2m], the energy dispersion of 
non-interacting free fermions with mass m), c^.^ and Ck;cr 
for the canonical creation and annihilation operators re- 
spectively for fermions with spin index a, gw{\\q\\) = 
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f)(||<7||) for the Fourier transform of the two-body in- 
teraction potential v{r — r') (assumed to be isotropic), 
and = 1/** for the volume of the (macroscopic) d- 
dimensional hyper-cubic box occupied by the system. 
The wave- vector sums in Eqs. (2) and (3) are over a reg- 
ular lattice (the underlying lattice constant being equal 
to 271 / L) covering in principle the entire reciprocal space. 
The Hamiltonian in Eq. (1) can be fruitfully employed to 
investigate the consequences of interaction and interplay 
between this and anisotropy (as implied by e^) on the 
observable quantities in many-particle systems. 
On effecting the replacements 



(4) 



and restricting the wave-vector sums to Nj^ wavevcc- 
tors uniformly distributed over the first Brillouin zone 
(IBZ) associated with the Bravais lattice spanned by 
{-RjIIj = l,...,Ni,}, the Hamiltonian in Eq. (1) goes 
over into the conventional single-band Hubbard Hamil- 
tonian [7-10] H corresponding to the on-site interaction 
energy U: in cases where k + q and p — q on the right- 
hand side (RHS) of Eq. (3) lie outside the IBZ, these 
vectors are to be identified with the vectors inside the 
IBZ obtained from k + q and p q through Umklapp 
processes. On relaxing the replacements in Eq. (4), while 
maintaining the above-mentioned restrictions concerning 
the wave vectors, the Hamiltonian in Eq. (1) coincides 
with an 'extended' Hubbard Hamiltonian. 

Along the lines considered in [6], it can be shown that 
for the Fermi surface Sp-a, pertaining to the uniform 
metallic GS of H, one has 



{kUe.-e<Jie+-e>J = 0}, 



where £+ = + 0+ and 



i4{k) 



(5) 



(6) 



On the basis of similar arguments as presented in [6], it 
can be shown that (sp — ^k-a) ~ ^ implies that (e+ — 
^fc-o-) ~ and vice versa so that, for fe e «SF;cr, e^.^ 
and £^.^ are up to an infinitesimal difference equal. In 
Eq. (6), 



,_ / n^(fe). 



n^{k), 



(7) 



I At At 



x(*JV;o|c^;^C^,_,_q,.o.,Cfe+q/;<,Cp/;^/|*JV;o), 



(8) 



(9) 



where |^jv;o) stands for the A/'-particle GS of il, 

na{ky. = {^N;0\Ck;aCk;a\^N;0) (10) 

for the GS momentum distribution function, and 

t?fc;a:=-(*iV;o|[4;a.[Hl,Cfc;,]-]^|*iV;o), (H) 

in which [ , ]_ and [ , ]+ denote commutation and anti- 
commutation respectively. Making use of the canonical 
anti-commutation relations for c^.^ and Ck;a, one readily 
obtains 

i?fe;. = ^tZ)(0)-l^zZ)(||p||)n.(fe + p) 



(for bounded w)(0)). (12) 



Note that N/il=:no is the total concentration of the 
fermions in the A^-particle GS of H. For the Hubbard 
Hamiltonian, making use of the prescription in Eq. (4) 
and the fact that J2p ncr(fc + p) = No- (recall that in the 
case of the Hubbard or an extended Hubbard Hamilto- 
nian the latter sum is restricted to points inside the IBZ), 
one obtains [6] 

"^kia = n — Ha = ria (Hubbard Hamiltonian), (13) 



where 



(14) 



stands for the number of fermions with spin index a per 
site in the uniform TV-particle GS of the Hubbard Hamil- 
tonian; thus n:=N/Ni^. In Eq. (13), a denotes the spin 
index complementary to a, that is for <j =|, a =[ and 
vice versa. 

In cases where i(}(||<7l|) is unbounded for ||q|| = 
(such as is the case for the two-body Coulomb poten- 
tial), the expression in Eq. (12) has to be regularized; 
this is achieved by formally taking into account the inter- 
action of the fermions with a positively charged uniform 
background through supplementing the Hamiltonian in 
Eq. (1) with gH', where 



H':=-now{K)N, k i 0, 

in which N stands for the total number operator, 
lowing this, one obtains 



(15) 



Fol- 



p 

(for unbounded w(0)), (16) 

where leaves out p = 0. For systems in the thermo- 
dynamic limit, sums of the form f^^^ X)q (o^ Sg) 
can be replaced by {2it)~'^ J A'^q. In any of the above 
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cases, whether one deals with the Hubbard Hamihonian 
or with a more general Hamiltonian, from the expressions 
in Eqs. (12) and (16) it is easily verified that 



(17) 



where T,^^{k) stands for the Hartree-Fock self-energy. 
From Eqs. (6), (7), (9) and (17), one readily deduces 
that (c/. Eq. (56) in [6]) 

^Ak)et,a + (1 - n.(fe))£>, = £fe + nS-(fe). (18) 

This result will prove useful in our later considerations 
in this paper. 

By introducing the decomposition 



(19) 



which implies no restriction so long as ncr{k) ^ 0, one 
can express e^.^, introduced in Eq. (6) above, as follows: 



(20) 



From this, Eq. (6) and the results in Eqs. (18) and (17), 
one similarly has 



^k;.=^k + g 



(21) 



From the expression in Eq. (5) above and the ensuing 
details (see the text subsequent to Eq. (6)) one infers the 
following relationships from the expressions in Eqs. (20) 
and (21) (for kp-c- S iSf;^) 



g 



(22) 

(23) 



In other words, to the leading order in (fe — feF;<7)) one 
has 



^k;a ^ ■^fe;<7 for fc ^ ^f;^- 



(24) 



Note that here, as well as in Eqs. (22) and (23), Sp-^a 
is considered as given, so that the expressions in these 
equations reflect some properties of Sp-cr rather than nec- 
essarily fully defining it. 

In analogy with the expression in Eq. (19) above, we 
introduce the atixiliary function rik-a, satisfying 



n<7(fe) Vk;a-='&k;^ (^F - Sfep, 



(25) 



This expression is general for all fe for which r\a-{k) ^ 0. 
In the light of Eq. (23), and since, for interacting GSs, 
ncr(fe) ^ in a neighbourhood of <Sf;<7 (see later), we must 
have 



%;o- ~ for fe ^ ^Fja- 

Defining (k;a as follows 

Ck;a-=^k;a (^F - £fcp;<,), 

g 

on account of Eq. (22) we similarly have 

Ck;a- ~ for fe — > 5f;<t- 



(26) 



(27) 



(28) 



On the basis of the above expressions, from Eqs. (20) and 
(21) for fe feFjCT e <Sf;<7 we deduce 

^k;a ~ M + 0'{kF;a) ' {k - k^;a) + g Ck;a, (29) 



^fe;a ~ /« + a{K;a) ■ {k - -gK{k) (Cfe;a - Vk;a), 

(30) 

(31) 



where 



a{kF;a)-= ^kek\k=k. 



Aa{k):= 



l-n<,(fe)' 



(32) 



In Eqs. (29) and (30), fj, stands for the 'chemical poten- 
tial' which for metallic states is infinitesimally greater 
than the Fermi energy £f (for details see [6,78]). 

Making use of the equation of motion for the operator 
Cfe;CT in the Heisenberg picture, that is 



d ^ 

iTl^Ck-ait) = [ck-a{t),H]_, 



(33) 



from the defining expression in Eq. (8), one obtains 

/3fc;a = -|^4Ga(fe;i-i') -£fen,(fe)}, (34) 
g '-''^ t'\.t 

where the single-particle Green function in the time do- 
main is defined as follows: 

G,{k;t- i'):= - i (*iV;o|T{cfe;<,(f)c^(t')}|*iV;o), (35) 

in which T stands for the fermion time-ordering operator. 
We have (below r? J. 0) 



h^G,{k;t-t') 



1 f-^ Ac 

1 

= -J deeA„{k;e), (36) 



where Acr(k;e) stands for the single-particle spectral 
function, defined in terms of Ga{k;z), the analytic con- 
tinuation of Ga{k;s) = lininio G a {k; s ± irj), £</U, as 
follows: 
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Aa{k; ^y-=^^{GAk; e-iri)- G^k; e + iri)], r? i 0. 

(37) 

With (see Eq. (10) above) 
1 f 

na{k) = T d£A„{k;e), (38) 

J -oo 

from Eqs. (19), (34), (36) and (38) we deduce that 



SY-y.={k II Efe + TiT,„{k; e^) = £f}, 



(45) 



}, (39) 



which in combination with Eq. (22) leads to 



j!l^deA,{k;e) 



Ef for k —^ kp-cr & Sj,-a ■ (40) 



With 

^fe,.:=na(fe-J-na(fe^J, (41) 
from Eqs. (38) and (40) we obtain 

^fe,. = ^ ^ [Mk^;a, £) - Mkt-y, e)] , (42) 

which is alternative to the standard expression 

Zk,._„ = ^ ds [A,(fe-,; £) - A„{k+^; s)] , (43) 

directly obtained from the spectral representation of 
r\a{k) in Eq. (38) above. The expression in Eq. (42) 
has been exphcitly derived in [6] for the uniform metalhc 

GSs of the single-band Hubbard Hamiltonian. Here as in 

> 

[6] the single-particle energies s^.^ (associated with some 
fictitious particles [6]) correspond to the single-particle 
spectral function Aa{k;e), for which we have [6] 

A.{k;s) = h{u<(k)6{e - e<.J + v>{k)6{e - £> J}. 

(44) 

It is readily verified that the exact Zk^_^ is obtained by re- 
placing Aa{k] c) in Eqs. (42) and (43) by Acr{k] e). In this 

> 

connection note that, for stable (uniform) GSs, e^.^ < n, 
Vfc [6]. The GS total energy of the interacting system is 
similarly exactly reproduced by replacing A(,{k] e) in the 
pertinent expression (see Eq. (51) in [6]) by A„{k;e). 



§ 3. EQUATIONS DETERMINING THE FERMI 
SURFACES OF UNIFORM METALLIC GROUND 
STATES 

The Fermi surface iSpjcr is conventionally defined as fol- 
lows [79,80,6] 



where, on general grounds, lm[S(j(A;; e^)] = 0, Vfc 
[81,6,78]. From Eqs. (17) and (23) we deduce 



£fe -I- nT^{k) = Ef for fe e <Sf;^. 



(46) 



Comparison of this expression with the defining expres- 
sion for iSf:ct in Eq. (45) reveals that to the leading order 
in (A; — k-p^a) we must have 

S<,(fe;eF)~S^"(fe) as fe ^ fcF;^ e ^f;^. (47) 

The significance of this expression becomes more ap- 
parent by considering the Kramers-Kronig relation for 
Re[S^(fe;e)] in terms of Im[S<^(fe; e)]; since [82,78] 

S,(fe;e)-E^^(fe) = o(l) for |£| ^ oo, (48) 

the above-mentioned Kramers-Kronig relation reads 



Re[S<,(fe;£)] =Sf (fc) 

de' sgn(/i - e') lm^^{k\ e') 



f°° de' 
-P — 

J-oo 



s' - e 



(49) 



From the leading term in the asymptotic series of the 
integral on the RHS of Eq. (49) for e — > Sf, making use 
of the fact that Im[I]o-(fc;£)] ~ as £ — > Sf, we obtain 



5],(fc;£F) = S«^(fe) + 5,(fe), 



(50) 



where 



1 (Ji£ 

S^{k):=- / — Im[E<,(fe; + s) + S<,(fe; £f - £)]. 

Jo £ 

(51) 

Comparison of Eq. (47) with Eq. (50) reveals the neces- 
sity of the condition 



S^{k) = 



(52) 



in order for k £ S^-^. In the light of the results in 
Eqs. (47) and (50), we thus arrive at the following alter- 
native defining expression for S^-^ {cf. Eq. (45) above): 

5f;^ = {fe II £fc + ?iS«^(fe) = £f A S^ik) = 0}. (53) 

The necessity of the additional condition So- (fe) = in the 
defining expression for S^-^a stems from the fact that the 
equation in Eq. (46) applies for all fe on a pre-supposed 
5f;<t so that, in general Eq. (46) cannot be sufficient fully 
to define S^-a ■ The significance of this aspect can be eas- 
iest appreciated by specializing Eq. (46) to the uniform 
metallic GSs of the single-band Hubbard Hamiltonian, 
in which case 'S^^{k) is independent of fe and equal to 
i[/nff (see Eqs. (12), (17) and (4) above). The suffi- 
ciency of £fc -|- ?iY,^^{k) — £f to define S^-a would in the 
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present case incorrectly imply that all solutions of the lat- 
ter equation, which constitute the entire si^^ (the Fermi 
surface corresponding to Sk and the partial number A^^- 
of fermions with spin index a in the exact interacting 
GS), were in general points located on Sj,-^^', the uncondi- 
tional equality S^-a = <sj,^aj rather than the relationship 
Sf-o- C Sp^l deduced in [6], would rule out the feasibil- 
ity that the pseudogap phenomenon is described in the 
space spanned by the uniform metallic GSs of the single- 
band Hubbard Hamiltonian [6], contradicting the results 
in [83-87]. Note that for isotropic GSs the condition 
Sa{k) = is automatically satisfied at any k for which 
Eq. (46) is satisfied, on account of the fact that the metal- 
lic nature of a GS implies that Sa-{ik) = for at least 
one k (and all points related to this by symmetry) on 
<Sf;(7, which by isotropy must remain valid for all k on 

We note that along the lines of reasoning as in [6] it 
can be shown that n(^(fe) is continuous (see Eq. (67) in 
[6]) at all k which satisfy Eq. (46) but fail to satisfy 
Eq. (52). We further note that, since for k G S^-a up 
to infinitesimal corrections [6] we have e^.^ = e^.^ = £f, 
Eq. (18) is seen to conform with the result in Eq. (46) 
above. 

It is interesting to rewrite the condition S„{k) = for 
fe e iSfjct in the following appealing form: 

Im[E.(fc;.)] ^ f-^^ Imp^^ ^ ^ ^^^^ 

(54) 

The integrands on both sides of this expression are non- 
negative over the pertinent ranges of integration. In the 
light of this and of the inherent symmetry of these inte- 
grands for both |e — EfI ^0 and |e — EfI ^ oo (to the 
leading orders in {e — Sp) and !/(£ — £f) respectively [78]), 
Eq. (54) can be viewed as implying some degree of simi- 
larity between |Im[Eo-(fc; £f — £)][ and |Im[Eo-(fc; £f + £)]| 
at intermediate values of £ (as measured with respect to 
£f) when fe e iSfjct; for fe ^ iSfj^, this similarity should 
be, if not absent, less pronounced. We point out that 



We believe that the apparent violation of the result in 
Eq. (54) by the first-order self-energy in terms of the screened 
interaction function (as evident from the middle figures in 
Figs. 16 and 17 in [88], corresponding to fc = l.Ok^) signifies 
the shortcomings of this approximate self-energy, as explicitly 
discussed in [78] (for instance, compare the first-order result 
in Eq. (278) with the exact result in Eq. (209) and in the 
light of the contrast between the two results, consider the 
expression concerning Im[So-(fe;£)] in Eq. (234a); in this con- 
nection, recall that in the paramagnetic state, dealt with here, 
no;9 — no;CT as encountered in the latter equation is equal to 
zero). 



since Scr{k) = applies for all fe e S^-^cr, in principle one 

has (sec later) 

Vfe'5,(fe')|fe,=feT =A,(fe^)n(fe) for fe e ^f;^, (55) 

where n(fe) stands for the (outward) unit vector normal 
to Sp-a at fe G Sp-a', here fe~ and fe^ denote vectors in- 
finitesimally close to fe, with the former inside and the 
latter outside the Fermi sea, and Acr(fe) stands for a scalar 
function which cannot be identically vanishing. For con- 
ciseness, in this work we shall denote the Fermi sea cor- 
responding to fermions with spin index a by FS<j and its 
complementary part with respect to the entire available 
reciprocal space, by FS^ [6]. 

In writing the expression in Eq. (55) we have taken 
account of the fact that Sa{k) may not be continuously 
differentiable in a neighbourhood of some or all points 
fc £ Sf-ct (whence fc^ in Xrr{k^)); however, in Eq. (55) 
we have not taken account of the possibility that Scr{k) 
may not be differentiable in a neighbourhood of fe e Sp-a. 
The latter possibility in fact arises when, for instance, 
fermions in d spatial dimensions interact through the 
long-range Coulomb potential, in which case VfcS"''(fe) 
is logarithmically divergent for k approaching any point 
feF;<7 on <Sf;<7 (both from inside and from outside the 
Fermi sea). It follows that, for metallic states of sys- 
tems of Coulomb-interacting fermions for which quasi- 
particles are well-defined in the neighbourhood of Sp-a, 
^ the above-mentioned singularity of VfcE™(fe) must be 
fully cancelled by a singular contribution arising from 
VfcS'<j(fc) (see Eq. (50) above), so that, for these metallic 
states, Eq. (55) cannot be meaningful. 

§ 4. EXPOSING A FURTHER CLOSE 
RELATIONSHIP BETWEEN Sa(fe;eF) AND Sf;^ 

The fact that, for fe ->■ k^.^ G S^-a, ^^(fe; £f) ~ ^^^''(fe) 
to leading order in (fc — fcF:cr) (see Eq. (47) above), de- 
serves some closer inspection. From Eqs. (12), (16) and 
(17) it is evident that S"''(fe) is fully determined by 
no-(fc). For uniform GSs this function is the Fourier trans- 
form of the single-particle density matrix ga{v — v'), or 

^<^(^) = j Wf"''^'^"'^^^- ^^^^ 

The singular nature of n^ik) for fe S iSfjo- has far-reaching 
consequences for the behaviour of Qaiv) for jjr]] — > oo 
[90,6] and consequently for that of S™(fe) for fe in a 



* This implying the differentiability of Ea(fe;ep) and con- 
tinuity of V)jSCT(fe;eF) in the infinitesimal neighbourhood of 
[6,82]. 
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neighbourhood of Sp-a- To illustrate this aspect, it is 
convenient to introduce 



= — where r:= 
r 



(57) 



Without loss of generality, we explicitly consider the case 
when d ~ 2 and assume that the angle ip in the cylindrical 
coordinates {k, ip) of k is measured with respect to the 
direction specified by 6^. We thus re- write Eq. (56) as 
follows: 



1 



where 



I4k,r):= d<^e''=''™«(^)n<.(fc,<^), 
Jo 



dkkl{k,r) 



(58) 



(59) 



in which ncr{k,(f) = na{k). Making use the stationary- 
phase method [89], we readily obtain 



laik, r) ~ 2\ — r\„(ker) cos(kr — 7r/4), r 
kr 



00, (60) 



where we have made use of the fact that by the time- 
reversal symmetry of the GS, 



no-(A;,7r) = n„{k,Q) = na{ker). 
Thus, from Eqs. (58) and (60) we have 
2 



(61) 



(27r)3/2ri/2 



dk na{ker) 
X cos{kr — TT / A) , r 



oo. (62) 

Let kp-a point in the direction of r, i.e. kp-^ = k^-^er. 
The singular nature of na{ker) at fc = k^i-a implies that 
for the further simplification of the expression on the 
RHS of Eq. (62) it is necessary [78] to expresses the k 
integral over [0, oo) in terms of at least two k integrals 
over the disjoint intervals [0, k^-a) and {k^-^cr, oo); this sub- 
division should be carried out for all k points at which 
n<j(fcer) is singular; however, to avoid any complication 
unnecessary to our present considerations, here we deal 
with the case where the most dominant singularity of 
r\a-{ker) is located at fc = kp-^a- Assuming r\cr{ker) to 
be discontinuous at fc = kp-,^ and continuous elsewhere, 
through integration by parts, from Eq. (62) we obtain 
[90,78] 



1 /2 

Zk^.^e^kp.g sin(fcF;o-r - 7r/4) 

2V27r3/2 



r3/2 



1 /2 

r\a{kzer)kz sm{kzr - tt / 4) 
2i/27r3/2 ;372 ' 



+ 



00, 



(63) 



where for systems defined on lattice k^er denotes the vec- 
tor from the origin to the Brillouin zone boundary (note 



that for these systems our assumption with regard to the 
uniformity of the underlying GSs implies that we should 
only consider r € {Rj\\j = 1,...,A^l}) and for those 
defined on the contimmm fc^ = oo so that nCT(fcze,.) = 0; 
for the former systems and for a given e^, n(j(fc2^r) may 
or may not be vanishing. In what follows, for simplic- 
ity we set r\a{k.^er) equal to zero. One observes that the 
full information concerning the shape and dimensions of 
<Sf;(7 is already contained in the leading-order term in the 
large- 11 r]] asymptotic series of ga{r). 

For the Fock part of the self-energy in the coordi- 
nate representation and for arbitrary GSs and two-body 
interaction functions v{r — r') we have [78] 



^^('^,■^0 = -^v{r -r')gcr{r',r), 



(64) 



where in the case of uniform GSs Qair' ,r) is to be re- 
placed by Qair' — r). It follows that the information 
concerning S-p-a contained in Qaifif') is similarly con- 
tained in ^ S^(r,r'). In particular, for any two-body 
potential other than those for which v{r — r') = 
when llr — r'|| > i? > 0, the behaviour of gu{r^r') for 
llr — r'll — > oo directly determines that of S^(r,r') for 
llr — r'll — > 00. For instance, for the uniform GSs of sys- 
tems of particles in two-dimensional space (i.e. for d = 2) 
interacting through the Coulomb potential and for which 
r\„{ker) is discontinuous at fc = kp-a, from Eqs. (63) and 
(64) wc have (following our above convention, for uni- 
form GSs we use the notation S^(r — r') for what in the 
general case is denoted by S^(r,r')) 



,1/2 



h (27r)3/2eo 



sin(fcji.;CTr - 7r/4) 

^5/2 ' 



00, 



(65) 



where stands for the particle-charge squared and eo 
for the vacuum permittivity (note that for — e the elec- 
tron charge, e'' jh = 3.874 . . . x 10^^ S is the quantized 
Hall conductance). ^'^ The fact that the oscillatory be- 



^ This observation corroborates the validity of the alterna- 
tive definition for iSf;ct in Eq. (53) above. 

^" We note that the direct evaluation of the expression in 
Eq. (65) from the Fourier integral representation of E^(r — r') 
is not straightforward. This is partly because Yf^(khr) is non- 
analytically singular at fc = fepio- (recall that in determining 
the expression in Eq. (63) we have explicitly assumed that 
the leading singularity of no-(fce,.) consists of a discontinuity 
at fc = fepio-) so that, for the determination of the large-l|r — r'lj 
asymptotic expansion of S^(r — r'), one has to employ the 
Laplace method [89] which requires analytic continuation of 
E^(fcer) into the complex fc plane. In this connection we 
point out that the appearance of sin(fcF;crr — 7r/4) on the RHS 
of Eq. (65), rather than cos{kF;<rr—'iT/4), should be a reminder 
of the non-triviality of the result in Eq. (65). 
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haviour of the leading term in the asymptotic series of 
T,^{rer) for r — > 00 is entirely determined by kp-^ in- 
dicates that, similar to r\„{k), T,^^{k) must be singular 
for all k on Sf;^', in the present case, where E^(rej.) 
decays with one power of 1/r faster than g(^{rer) for 
r ^ 00, the assumed discontinuity of n^{k) at fe = kp-a^r 
implies singiilarity (explicitly, a logarithmic divergence 
in the behaviour of VfcE^(A;) at fc = kp-cr^r- These 
facts expose a fundamental distinction between systems 
of fermions interacting through potentials of finite range 
and potentials of infinite range. Notably, whereas for the 
uniform metallic GSs of the single-band Hubbard Hamil- 
tonian one has n„{k) > i for k inside the Fermi sea and 
infinitesimally close to iSf;^ [6], as we show later in this 
paper (§ 6), this is not necessarily the case when the inter- 
action potential is of longer range than the intra-atomic 
potential encountered in the Hubbard Hamiltonian. In 
this connection, recall that for the uniform metallic GSs 
of the conventional single-band Hubbard Hamiltonian, 
E^''(A;) = j-Un^, which is independent of k and thus 
regular for all k £ IBZ. 

It is important to point out that n(j(fe) is singular {not 
necessarily discontinuous) at all k satisfying Eq. (46) 
above. However, since according to Eq. (53) the sat- 
isfaction of Eq. (46) at a given k is not sufficient for 
fe £ Sf-o-, it follows that the fe that for a given direction 
determines the leading term in the large-r asymptotic se- 
ries of Qa{rer) (and of E^(rer)) is not necessarily a point 
of iSf;<t. Since such a fe point satisfies Eq. (46) but fails 
to satisfy Eq. (52), it follows that at this point 'n^ik) is 
continuous (see the paragraph preceding that containing 
Eq. (54) above) whereby the leading term in the large-r 
asymptotic series of Q^{rer) (and of E^(rer)) acquires 
an anomalous exponent, leading to a stronger power-law 
decay in the mentioned asymptotic term in comparison 
with the case corresponding to Zke^ ^ for ker G S^-a', 
with reference to the expressions in Eqs. (63) and (65) 
above (where k = kp-^), the above- indicated anomalous 
exponent would read 3/2-1-7 and 5/2-1-7 respectively, 
with 7 > 0. 



This statement can also be explicitly verified from the 
expression for S"''(fc) = {'g/K)'dk;<T (see Eq. (17)), making use 
of Eq. (16). 

For some pertinent details we refer the reader to an ex- 
tensive analysis of the largc-r asymptotic behaviour of Qair) 
pertaining to uniform and isotropic GSs in appendix J of [78]. 
From this analysis and the bounded result in Eq. (65) above 
one can readily infer the nature of the singularity (i.e. loga- 
rithmic divergence) of VfcE^(fc) at fe € Sr-a- 

^'^ The analysis underlying this statement is analogous to 
that presented in [6] for the case of the single-band Hubbard 
Hamiltonian (see Eq. (68) and the subsequent text in [6]). 



§ 5. GROUND-STATE EXCHANGE AND 
CORRELATION POTENTIALS: A TEST CASE 

A significant aspect of the result in Eq. (46) is made 
explicit as follows. For concreteness, here we consider the 
uniform GSs of systems defined on the continuum; these 
uniform GSs are necessarily isotropic. Since lattice mod- 
els are excluded from our considerations in this Section 
(in a future work we explicitly deal with the more deli- 
cate case of lattice models), here we relax the definition 
for Ua- introduced in Eq. (14) above and introduce 



(66) 



Ua, the total concentration of the fermions 



With n:=n. 

in the A/'-particle uniform GS of the system, we define the 
spin magnetization fraction m as follows 



ria - ria 



m:=- 



(67) 



Denoting by E{N, M) the GS total energy of the system 
described by the H in Eq. (1), corresponding to the total 
number of particles N = fin and the total spin magneti- 
zation M = Q{ncr — ria), we define 



£{n,m):=^E{N,M), 



(68) 



which we assume to be a sufficiently smooth function of 
n and m in the thermodynamic limit [91] (appendix B 
herein). In this limit, for a metallic GS we have 



d_ 
dn 



(n£{n, to)) . 



(69) 



From this we obtain 



d . > d£{n,m) 

{n£{n,m))=eF-{mTl)^^, a 



We now decompose £{n,m) as follows: 

£{n, to) = £k(n, m) + fxc(?^, rn), 



(70) 



(71) 



where £\^{n,m), the 'kinetic energy' contribution to 

£{n,m), corresponds to the GS of Hq in Eq. (80) be- 
low whose associated n and to are constrained to be the 
same as those pertaining to the exact GS of H. Thus, 
Eq. (71) defines the exchange-correlation energy contri- 
bution £:xic{n,m). By the assumption of the isotropy of 
the GSs under consideration and in consequence of a Lut- 
tinger theorem [92,80], according to which the number of 
fe points enclosed by Sp-a- and S^^^ are equal, we have 
Sp-cr = Sjr^a^ Vcr. Prom this and in analogy with the re- 
sult in Eq. (70) we have 
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d 



[n£i,{n,m)) =Sk,,^-{mTl) , a=|^ 



(72) 



Combining the expressions in Eqs. (70), (71) and (72), for 
the exchange-correlation potential ^i^c-a, defined through 



IJ-^c;a~-^{ri£^c{n,m)), ae{a,a}, (73) 
we deduce that 

/«xc;a = £f - £fep;„ - (tot 1) , " = \ ^ V^) 

By subtracting y^xcis from /ixc;(7, employing the expres- 
sion in Eq. (74), we obtain 



d£^c{n,m) _ 1 
dm ~ 2 



where 



efe;a:=£fe + Mxc;a, " € {ct, a). 



(75) 



(76) 



Prom Eq. (75) we infer that in a paramagnetic state, 
where ek-„ = ek;s, we must have 



dm 



= (for paramagnetic states). (77) 



Prom Eqs. (74) and (77), for paramagnetic states we thus 
obtain 

(78) 

which is a well-known theorem due to Seitz [93]. 

A result which from the perspective of our present work 
is of special interest is 

Mxc;a = ?iS-(fcF;a) " (to T 1)^^^, a = { 

(79) 

which is deduced through comparing Eqs. (46) and (74). 
Por completeness, the non-interacting Hamiltonian 



(80) 



with €k;a as defined in Eq. (76) above, is the so-called 
Kohn-Sham Hamiltonian [76] encountered within the 
framework of the GS density-functional theory due to 
Hohenberg and Kohn [75] (for a review see [77]). 

By repeating the above considerations in terms of the 
GS energy per particle £{n) (an explicit function of solely 
the total concentration n and not of the spin magnetiza- 
tion fraction m), and thus by introducing, analogously 



as in Eq. (71) above, the energy functions £y:{n) and 
£xc{n) satisfying £i^{n) +£xc{n) = £{n), we arrive at (c/. 
Eq. (78)) 



-(0) 



(81) 



in which Cp"^ stands for the Permi energy of the under- 
lying paramagnetic non-interacting A^-particle GS, and 
(c/. Eq. (73)) 



(82) 



It is observed that, for paramagnetic GSs, Eq. (81) co- 
incides with Eq. (78). In the case at hand, the non- 
interacting Kohn-Sham Hamiltonian Ho in Eq. (80) is 
defined in terms of the spin-degenerate energy dispersion 
(c/. Eq. (76) above) 



(83) 



where /Uxc is defined in Eq. (82). In the light of Eq. (81), 
the result in Eq. (83) expresses the fact that the energy 
eigenvalue corresponding to the highest occupied single- 
particle Kohn-Sham state coincides with the exact chem- 
ical potential. Following Eq. (77), in the present case the 
counterpart of the expression in Eq. (79) is as follows: 



Mx 



(84) 



Note that, since /Xxc (A*xc;cr) is a potential, it is determined 
up to an additive constant, explicitly independent of n 
(n and m); the difference /Xxc;<7 — Mxc;a- is, however, not 
arbitrary. 

We have employed the expression in Eq. (84) special- 
ized to the paramagnetic GS of the Coulomb-interacting 
uniform-electron-gas system for = 3, with Sk = 
?i^||fe|p/[2TOc], where nic stands for the bare electron 
mass. Decomposing /Uxc;o- into exchange and correlation 
contributions, that is 



/^xc;(T — /^x;cr T /^c;(T, 



(85) 



in Hartree atomic units, where g is equal to unity, we 
have 



/ix-<7 = / dx x\n 



1-x 



Mc;<7 



70 



' s Jo 



dx xln 



1 + x 
l~x 



n<-^\hx), Va, (86) 



1 + x 

x[r^^{hx)-n^'^ihx)], Va; (87) 



Here and in the subsequent expressions wo use the spin 
index a in order to avoid the confusion that often ensues 
on suppression of a; some workers use n(fe) to denote r\a{k) 
for paramagnetic GSs, while others employ n(fe) to denote 

E."4fc)- 
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here the bars (such as in k^) indicate that the respec- 
tive quantities arc in Hartree atomic units. In Eqs. (86) 
and (87), n'^\kpx) stands for the momentum distribu- 
tion function pertaining to the GS of the non-interacting 
electron-gas system (equal to 1 for < x < 1 and for 
X > 1) and Ts for the average inter-particle distance ro in 
units of the Bohr radius; 



7o:= 



97r 



1/3 



1.919... 



(88) 



and fcp = llfepjo-ll denotes the Fermi wave number in units 
of the inverse of the Bohr radius, i.e. 



rvp — 



Making use of the standard result 



Jo 



dx a; In 



1-x 



1 + x 



(89) 



(90) 



from the expression in Eq. (86) we trivially obtain the 
well-known exact result 



Mx;<T = 



-7o/7r -0.610887. 



We note that for d = 3 (see Eq. (82) above) 

Mxc;(T — ^xc "7^ \ ^xc* 

3 dr. 



(Hartree). (91) 



(92) 



With fxc = £^ + £c {cf. Eq. (85) above), in which [94] 
s -370 



47r Ts 



(93) 



from Eq. (92) we identically reproduce the result in 
Eq. (91). 

In order to display the significance of the expression 
for /Xc;cr in Eq. (87), we employ the ncr(fc) due to Daniel 
and Vosko [95] evaluated within the framework of the 
random-phase approximation (RPA). Our numerical 



Our numerical calculations reveal that for rs > 



n"'*(fcF;<.) < h (foi- rs = 6.09887. 



(fc^.) = 



n 



0.182397 . 



6.09887 . 

0.31760 ... so that for this rg we have Z^^ ^ 
and [n:^'"*(fc-,) + nf "(fc+<,)]/2 = 0.408801 !. .) and that for 
rs > 7.769269 . . ., nf-\ki.J) < nf''{k+.„). The latter of course 
implies that, according to the RPA, the paramagnetic GS 
of the uniform electron-gas system should be unstable for 
rs > 7.769269 .... For completeness, according to Bloch [96] 
(see [91], pp. 682-684), within the framework of the Hartree- 
Fock approximation the latter GS is unstable towards a uni- 
form ferromagnetic state for rs > 5.4502186 .... This is sug- 
gestive of the possibility that the instability as predicted by 
the RPA (which is an artefact of this approximation) may 
also be one towards a ferromagnetic GS. 



calculation of flc-.a based on the expression in Eq. (87) 
and the latter GS momentum distribution function re- 
produces, to within the numerical accuracy of our calcu- 
lation, the exact result [94] 

Mc:<T - — ^ In(rs) (Hartree) for r, j (94) 

(Figs. 1 and 2). This result implies that for | the 
integral on the RHS of Eq. (87) approaches zero like 
Ts In(rs). The fact that the latter integral should be van- 
ishing for Ts i can be understood by the observation 
that Tg is proportional to the coupling constant of inter- 
action so that, for rs i 0, n^{k) — n'^\k) should indeed 
approach zero for all k. However, the above-mentioned 
behaviour of the integral on the RHS of Eq. (87) very 
crucially depends on the precise way in which ncr{k) ap- 
proaches n'^^k) for rs J, 0; even the slightest deviation 
of r\„{k) from the exact n^{k) for rs i can signifi- 
cantly alter the behaviour of the latter integral for rs i 
and consequently of p,c;cr- We suggest therefore that any 
parametrized form of n„{k) should be required to repro- 
duce a pre-determined p,c;a through the expression in 
Eq. (87). To illustrate the significance of this aspect, by 
employing the expression in Eq. (92) above we have de- 
termined jj,c;<T through use of an interpolation expression 
for f c;cr [97,98] based on quantum Monte Carlo calcula- 
tions [99] for the Coulomb-interacting uniform-electron- 
gas system at a number of densities. Comparison of this 
j2c;a with the fi^-^r obtained through Eq. (87) in conjunc- 
tion with n^'''^(fc), reveals that to the numerical accuracy 
of our calculations the two results up to an unimpor- 
tant additive constant (independent of rs) coincide for 
small values of rg where n^''^(fc) — n'^\k) constitutes the 
leading asymptotic contribution to nc^{k) — n'^\k) (see 
Fig. 1). We have also determined jlc;cr by employing a 
very recent parametrized expression for ncr(fc), asserted 
by Gori-Giori and Ziesche [100] to be valid in the range 
Ts ^ 12, and observed no resemblance whatever between 
this /ic;cr and the aforementioned quantum-Monte-Carlo- 
based j2c;(T over the entire range rg ^ 12 (see caption of 
Fig. 2). The reason for this lies in the fact that, even at 
rs = 1, fic:a is smaller by one order of magnitude than 
Mx;<7 5 with the disparity between these two quantities fur- 
ther increasing for smaller values of rs (this as evidenced 
by the fact that fic-.a / P-x-.a for | 0); it follows that 
no parameterization of ncr(fc) is likely to reproduce j2c;o- 
through the expression in Eq. (87) unless this reproduc- 
tion has explicitly been enforced in the determination of 
the parameters of n„{k). 
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§ 6. BEHAVIOUR OF n<,(fe) FOR k 
INFINITESIMALLY CLOSE TO Sp^a 

Wc now proceed with the investigation of the be- 
haviour of r\c{k) for k approaching S^-a- Along the hnes 
presented in [6], for k G Sp-a we have 



A,(fe) = r,(fe), 



(95) 



where Ao-(fe) has been defined in Eq. (32) above and 

r,(fe):=^^^. (96) 



'fe;o- 



Making use of the expressions in Eqs. (29) and (30) 
above, we deduce that 



r<.(fe) 



-a(fcF;a) ■ (fc - k^-a) -gCk-o- 

a{kp-a) ■ {k - K-a) - g A^(fe)(C/e;a " Vk;a) ' 

A; ^ S^.,a, (97) 



where on account of the assumed stabiUty of the GS of 
the system both the numerator and the denominator on 
the RHS must be non-negative. The fundamental differ- 
ence between the expression on the RHS of Eq. (97) and 
its counterpart specific to the uniform metallic GSs of 
the conventional single-band Hubbard Hamiltonian (c/. 
Eq. (76) in [6]) arises from r]k;cr which for the latter GSs 
is identically vanishing. Since we have earlier considered 
the Hubbard Hamiltonian in some detail [6] , in what fol- 
lows, unless we explicitly indicate otherwise, ijkia ^ 0. 
As we shall see below, an identically non- vanishing rjk-a 
has far-reaching consequences for the behaviour of no-(fe) 
for k approaching S^-a- 

With reference to the expressions in Eqs. (28) and (26) 
above, we consider the following behaviours for (k;a and 

Vk;(T as k > ^F;(T ^ ^F\(t' 



\Ck;a\ - \Bik,,.)\ ||fe-fep.,|r, 

|Cfc;-| ~ \B{k^.,a) lnllfe-fep;,|| 



(98) 
(99) 



\r]k;a\~\C{k^;a)\\\k-k^.,4\ (100) 
\Vk;a\ ~ |C(fep;^) In life - fep;^||| life - (101) 

where B{k^-„), 7, C{kF-„) and r are parameters. Be- 
low we consider a variety of possibilities corresponding 
to various ranges of values for the parameters 7 and r. 
Since in the present case, in addition to the intervals in 
which 7 and r lie (0 < 7 < 1 or 7 > 1, for instance), the 
precise relationship between 7 and t (i.e. 7 < r, 7 = r 
and 7 > r) is significant, an exhaustive analysis of all 
possible cases will be too extensive to consider here. We 
shall therefore explicitly deal with a limited number of 
cases which arc representative and further pave the way 
for future complementary considerations. Below we con- 
sider the following cases: 0<7<t<1;0<t<7<1; 



0<r = 7<l;7 = l,T=l. For obvious reasons, when 
for instance Eq. (98) applies, the cases corresponding to 
< 7 < 1 and 7=1 should be considered separately 
whereas, when Eq. (99) applies, we need to deal with 
one case only, corresponding to < 7 < 1. Similarly for 
Eqs. (100) and (101) and < r < 1, r = 1 and < r < 1 
respectively. Note that in principle the values of B, C, 7 
and T can depend on whether fc approaches Sp-a- from the 
inside (outside) of the Fermi sea, in which case B{k^-cr), 
C{k^-a), 7 and r in Eqs. (98) - (101) should be denoted 
by B{k-,), C(fe-J, 7- and r" (i?(fc+J, C(fe+J, 7+ 
and T+) respectively. 



6.1. The case when 0<7<1, 0<t<1 



In this case from Eq. (97) we obtain 

1 Ck-c 



Taik) 



, k^S^,„. (102) 



In the present case, the assumption with regard to the 
stability of the GS of the system under consideration im- 
plies that 

Cfe;o- < and Cfe;o- < ??fe;(7, for fe — > S^-a, (103) 

must hold; violation of these conditions imply that for fe 
sufficiently close to S^-cr, would increase above ^, and 
^k <7 would decrease below respectively. 

To proceed, we need to consider three different cases, 
corresponding to 7 < t, to 7 > r and to 7 = t. In the 
last case, it is relevant whether Ck\a satisfies Eq. (98) or 
Eq. (99) and similarly whether rjk-„ satisfies Eq. (100) or 
Eq. (101). The case corresponding to 7 = r = 1, with 
Cfe;<7 and rjk-a satisfying Eqs. (99) and (101) respectively, 
concerns in particular (i.e. a priori not exclusively) sys- 
tems of Coulomb-interacting fermions (see the paragraph 
subsequent to that containing Eq. (55) above). 



6.1.1. 0<7<r<l 



For 7 < r, from Eq. (102) we obtain 
1 



ra(fc) 



fc > S^-^a^ 



K{k) 

which in conjunction with Eqs. (32) and (95) yields 



(104) 



(105) 



Note that in this case, where |Cfe:cr| > \rik:a\ for k suffi- 
ciently close to 5f;(7, satisfaction of the ffist requirement 
in Eq. (103) above implies satisfaction of the second re- 
quirement. With 
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na(feF;a):=2 [na(fe,. J + n,(fe+ J] , (106) 
from Eq. (105) wc obtain (see Eq. (41) above) 

Mk.;a) = l, ^fe,. =0 (7<t). (107) 

It follows that in the case under consideration the single- 
particle spectral function Aa-{kF-o-: e) is free from a quasi- 
particle peak at e = Ef (breakdown of the FL picture). 



6.1.2. 0<r<7<l 

For 7 > r, from Eq. (102) we have 

~ — +°°' k^S^.,^. (108) 

Since here |?7fc;cr| > |Cfc;cr| for ^ sufficiently close to S^-^a^ 
for such k the stability conditions in Eq. (103) imply that 
Vk-a > 0; hence we have +oo on the RHS of Eq. (108). 
The result in Eq. (108) leads to 



n^{kf.J = l, Zk,.^=0 (7>t). 



(109) 



In the case when such behaviour for n„{k) is unlikely, one 
is left to conclude that 0<T<7<lis equally unlikely 
to be satisfied for a realistic model of interacting fermion 
systems. It is interesting, however, to note that (see our 
pertinent discussions in § 1) n^{k^.^) = 1 is not dissimilar 
to that observed in the experimentally determined na-{k) 
corresponding to Bi2212 and Bi2201, at the Fcrmi-surface 
points adjacent to the M points of the Brillouin zone (see 
Figs. 2(gj and (i) in [73]). 



6.1.3. < r = 7 < 1 (includes Fermi liquids) 

Here we explicitly deal with the instances where the 
expressions in Eqs. (99) and (101) apply. As we have in- 
dicated earlier (see the paragraph following that includ- 
ing Eq. (55) above), for the metallic states of fermions 
interacting through the long-range Coulomb potential, 
VfcSf (fc) = {g/h)Vkdk;a (see Eq. (17)) is logarithmi- 
cally divergent for k = k^-^ € S^-„. It follows there- 
fore that the FL metallic states of Coulomb-interacting 
fermions fall into the category of systems dealt with in 
this Section. 



EYom the defining expression in Eq. (25) and the fact that 
%;ct ~ for k —f ki,-cr £ Sf-ct (the latter rendering an assump- 
tion with regard to the boundedness of 'Vkna{k) for k — > fcp;c7- 
redundant) it follows that a possible (logarithmic) divergence 
of VfciJfc.o. for k fcriCT directly implies an analogous diver- 
gence of 'VkVk a for k — > fepiCT- 



From the expressions in Eqs. (99) and (101) it follows 
that to the leading order in (fc — fepio-) we can write 



(110) 



where A~ and A+ stand for finite constants specific to 
fc = fepicr- Thus, in view of the expressions in Eqs. (99) 
and (101), for r„{k), as presented in Eq. (97) above, we 
have 



r.(fe) 



-gCfc;. 



1 



-{l-Xf)gAa{k)Ck;a (l-Aj)A,(fc)' 



(111) 



The requirement of the stability of the GS under consid- 
eration implies that 



A? < 1. 



(112) 



Solving Eq. (95) for fc = fc^^^, kp-a- G Sp-a, with r<^(fc) 
therein replaced by the asymptotic expression on the 
RHS of Eq. (Ill), we obtain 



Aa(fe,^J 



1 



which in conjunction with Eq. (32) yields 



(113) 



(114) 



1 + Vl - A? 
Prom this and Eq. (106) we have 

n-.(fc.;.) = ^( ^— + ^=t)' (115) 

2^1 + Vl- A^ 1 + Vl-A+^ 

which is not necessarily equal to 5. 

To gain insight into the quantitative values to be ex- 
pected from A+ and A~ . we recall that in the weak- 
coupling limit for the Coulomb-interacting uniform- 
electron-gas system (and thus in the paramagnetic phase 
of this system) one has ^"^ 



"<7(fcF;a)~l-"<7a 



for Ts i 0, 



where 



a:= 



7r^7o 



(116) 



(117) 



Here we are relying on the fact that for I the leading- 
order asymptotic behaviour of the exact nrr{k) — n^''(fc) is 
exactly reproduced by n^'''^(fe) — n^°'(fe), Vfe. 
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in which 70 has been defined in Eq. (88) above. The 
constants in Eq. (116) are in principle bounded func- 
tions of Ts; according to Daniel and Vosko [95], within 
the framework of the RPA, a~ = a+ w 1.7. As can 
be observed from the results in Fig. 3, are not con- 
stants; although for Ts approaching zero, a~ and a+ tend 
towards the same value (which we denote by a^), for 
sufficiently small rg, monotonically decreases towards 
zero, resulting in the fact that for Vs [ 0, no-(feF;<7)) as 
defined in Eq. (106) above, approaches \ more rapidly 
than asserted in [95]. 

Making use of the results in Eqs. (116) and (117), from 
Eq. (114) we obtain 



A- -1 



At ~1 



where 



a„a 



1 -Qqa 



1 - (a<xa)' 



\ (118) 



~ aa- for i 0. 



(119) 



We observe that in the regime of weak interaction, corre- 
sponding to fs -C 1, A~ t 1 and A+ [ —00. We point out 
that similar expressions to Eqs. (114), (116) and (118) 
apply to anisotropic metallic GSs, with the parameters 
therein dependent on the direction in the k space along 
which <Sf;ct is approached. 

From the above considerations it follows that in the 
weak-coupling hmit, corresponding to < a^a ^ 1, for 
fe — »■ >Sf;ct one must have (see Eqs. (22) and (23) above) 



1 + a^a, for k e FS^, 
i?fe;a-(eF-efe,;J/g \ -a„3, for fe e FS,,. 

(120) 



In arriving at this result we have made use of Eqs. (27), 
(25), (110), (116) and (118). 



6.2. The case when 7=1, t = 1: general 

Here we explicitly deal with the instances where the 
expressions in Eq. (98) and (100) apply. In the present 
case, for fe kp-a € Sp-cr we have ^® 



Cfe;a~ba(fe^J-(fe-feF;.), ^ ^ { (121) 



m;a Ca{kf.„) ■ {k - kp-a), ^ ^ {^'" 



With 



d^:=d^{kf.„) •n(feF;a), 

in which 

dAkf.^^):=Kikf.^) - c,{kf.,) 
from Eq. (97) we obtain 



(122) 



(123) 



(124) 



(125) 



r.(fe) 



k kp-a e ^f;^. (126) 



The assumed stability of the GS of the system under 
consideration implies that both the numerator and the 
denominator of the function on the RHS of Eq. (97) must 
be non-negative. From this, considering Oo- > 0, g > 
and A^{k) > 0, we arrive at the following conditions 



b->-^<o, 



>o. 



bt<--< 0, 
g 



>o. 



^ ^ gAt 



(I) 
(11) 
(III) 

(IV) 



(127) 



where 



K^:=K{k%a). 



(128) 



Making use of the expression in Eq. (126) above, the 
equation for Sy-^ in Eq- (95) reduces to the following 
quadratic equation for AJ: 



gd^{kjf- a, A J - (a, +g6T) = 0. 



(129) 



For interaction potentials with non-vanishing range, in 
general ^ 6J (see Eq. (125) above), whereby Eq. (129) 



Here as in [6], for k — » kr-^a € &;ct we assume (with- 
out loss of generality) tliat fc = k^-a T 11^ ~ kp-cr\\n{kF:a), 
k G FSct, FS(T, where n{kF;a) stands for the outward unit vec- 
tor normal to &;(t at kw-,^- Consequently, in what follows 
kf.„:=kF-a T Kn{kF-a), with k J. 0. See Eqs. (88) and (89) in 
[6]'. 



Although in prmciplc ao- = 0, g = and Aa-{k) = are 
feasible (separately or in combination), the cases correspond- 
ing to these possibilities are most conveniently dealt with by 
considering the limits of the pertinent results, to be presented 
in the following, for o<t i 0, g j and Aa{k) J, 0. 
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acquires a rich spectrum of solutions. Considering the 
second-order equation ax^ + hx + c = 0, for which we 
have the solutions x^ = —b/[2a] =F - 4ac)^/^/[2a], 
below in analogy we denote the possible two solutions 
of Eq. (129) by A^^. From Eq. (129) we immediately 
obtain the following results 



(130) 
(131) 



We note that for x± = u ± iv, with u and v real, we 
have x+x^ = + > from which it follows that 
in cases where the RHS of Eq. (131) is negative, the so- 
lutions (whether Aj^_|_ and Aj^_ or A^^,,, and Aj^_) can- 
not be complex. In connection with x± = u ± iv, we 
point out that, since g, and are real, the possi- 

ble complex solutions of Eq. (129) must indeed occur in 
complex-conjugate pairs. 

Below we separately consider A^ and A+. In our fol- 
lowing considerations we shall encounter which in 
analogy with 6J and d^ in Eqs. (123) and (124) respec- 
tively are defined according to 



cj: — Co-(fe^o-) ■ '^(^F;o-)- 



6.2.1. Considering na{k^.^) = A„ /(I + A^ ) 

According to (II) in Eq. (127) above, d^ > 0. Since 
according to (I) in Eq. (127), Ua + gb^ > 0, in view 
of d~ > from Eq. (131) it follows that A^_ A^_|_ < 
0. Thus, of the possible two solutions of Eq. (129) for 
A~ (which by our above argument, presented subsequent 
to Eq. (131), cannot be complex), one is positive and 
one is negative; with reference to Eq. (32) above, it is 
the positive solution, the 'physical' solution, that is of 
interest to our considerations that follow. Introducing 



C±{x,y)- 



1 

2^ 



"l±(l + 4a;(l + 2/))'/'], (133) 

for the physical solution of Eq. (129) for A~ we have 

A- =A-+=/:+{x-,y-), (134) 

where 

- gd- _ Kb-„ 



x„ :=- 



Va 



5x'„ 



in which (see Eq. (125) above) 



5x„:= 



gCa 



(135) 



(136) 



From (I) in Eq. (127) we deduce the following two iden- 
tical relationships 



Va > -1 



l + x-+5x-> 0. 



(137) 



From Eqs. (130) and (135) and the fact that A~ = A^_,_ 
(see Eq. (134)) and A~_ < 0, we further have 



x„ = 



K+ + K 



Aa 



Making use of the exact result 

£,(.,. + ^.) = l[l + 2(.+ l)(l + ^-|^ 



(138) 



1/2- 



(139) 



from Eq. (134) we deduce the following inequalities 

on x~ > 0) 



a; 51 



which after some manipulations give rise to 



a;^i 



2. 



(132) From Eqs. (136), (141) and (32) we thus obtain 



A-^1 



(141) 



(142) 



Since —2acr/g < 0, it follows that in the event c~ = (as 
is the case with the single-band Hubbard Hamiltonian 
^° [6]) we have A~ > 1 or ncr(fc~^) > \. The more 
general result in Eq. (142), allowing for r\a{k'^.^) < i, 
signifies the considerable influence that the (long) range 
of the interaction potential can have on the behaviour of 



For instance, at half- filling, for (i = 1 and £k = — 2tcos(fc), 
t > (here k is in units of the inverse lattice constant so 
that at half-filling ±/cF;a- = ±7r/2), in the limit of large U/t 
one has the exact result [101] (see also [102]) ncr{k) ^ ^ + 
4 ln(2) cos(fc)t/[/, k G [— tt.tt); it is seen that in this limit 
indeed no-(fc) > | for |fc| < fepio- and <na{k) < i for jfcj > kp-cr- 
Similar behaviours obtain for d = 2 (square lattice) and d = 3 
(simple cubic lattice) [101]. We point out that for d = 1 and 
Sk = — 2<;cos(fe), t > 0, following the exact Bethe-Ansatz 
solution due to Liob and Wu [103], at half-filling the GS is 
insulating for all U > 0, from which it follows that for the case 
at hand the absence of discontinuity in ncr{k) at fc = ±7r/2 
is not specific to U/t oo, but to all U/t > (in fact, 
in this case and so long as U/t > 0, d'"no-(fe)/dfe'" can be 
shown to be bounded at fc = ±7r/2 for any finite value of m). 
Conversely, contrary to the statement in [101], the absence of 
the usual discontinuity in na{k) (with reference to the above 
case, at fe = ±7r/2) does not necessarily imply absence of 
Fermi surface and thus an insulating underlying GS. 
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6.2.2. Considering n^{k+.„) = A+/(l + A+) 

In contrast with d~ which is positive, d+ has no a 
priori sign; we know only that (sec (IV) in Eq. (127) 
above) d+ < Oo-/(gA+). This necessitates us to consider 
separately the cases rf+ > and rf+ < 0. Here we have 
however (sec (III) in Eq. (127) above) -{a^ + gK) > 0- 
Thus following Eqs. (130) and (131), for g > we have 



A+ + A+ > n A+ A+ > n 



d+>0, (143) 



that is, the two solutions of Eq. (129) for A+ are posi- 
tive when d+ > while they have different signs when 
rf+ < (note that, by the argument presented following 
Eq. (131) above, in the case whore q!+ < the solu- 
tions cannot be complex); in the latter case, the nega- 
tive solution has a greater magnitude than the positive 
solution; in the former case (d+ > 0), where there are 
two positive solutions, we need to invoke the restriction 
A^ > A+ (see Eq. (153) below) in order to guarantee 
na{k~.^) — na{k^.^)=:Zk^^^^ > 0. In this connection, we 
note that for the case when d+ < 0, where we have 
the single positive solution A^_ for A+, the condition 
A~ > A+ is similarly not trivially fulfilled so that it nec- 
essarily restricts the ranges of variation of (some of) the 
parameters of the problem at hand. 



6.2.2.1. The case when d+ < 

Here, in view of the requirement Ao-(fc) > 0, A+ is to be 
identified with Ajr_ (following Eq. (143), in the present 
case AJ^ < 0), that is 

A+ = A+_^£_(a;+,2/+), (144) 

where C-{x, y) has been defined in Eq. (133) above, and 

xt:=^, yt:=^^xt + 5xt. (145) 



in which 



5xt:-- 



gc+ 



(146) 



Prom (III) in Eq. (127) we deduce the following identical 
relationships (c/. Eq. (137) above) 



< -1 ^ l + xt + 5xt < 0. 



(147) 



In addition, from Eqs. (130) and (145) and the fact that 
A+ = A+_ (see Eq. (144)) and A+^ < 0, we have 

(148) 



A+++A+_' ^-^Kf 
Making use of the exact result 



C-{x,x + 5x) = — 
2x 



1 - 2 



t5x \ 1/2 



(x+ 1/2)2 



(149) 



from Eq. (144) we arrive at the following inequalities 
(based on a;+ < 0): 



A+ > 1 



1 / T+r^T+ \ 1/2 



{xt + l/2fl < 



which after some manipulations give rise to 



A+>1 ^ 5xt< -2. 



From Eqs. (146), (151) and (32) we thus obtain 



(150) 



(151) 



c+ < ^ ^ A+ > 1 ^ n.(fe+J > \. (152) 

The fact that here the condition r\^{k'^.^) > ^ is in prin- 
ciple feasible clearly signifies the considerable influence 
that a non-vanishing range of the two-body interaction 
potential can have on the behaviour of n„{k'^.^). The 
requirement ncr{k^.^) — na{k'^.„)=:Zk^.^ > implies the 
necessity for the following condition to be satisfied: 



^^^±fe^>l for Sx-< 
At C-{x+,y+)- 



2 A Sx+> 



2. 



(153) 



In connection with the above, we mention that the com- 
bination Sx~ > —2, 5x+ > —2 corresponds to n^ikp.^) > 
^, r\„{kp.^) < i so that it does not imply any further re- 
striction; on the other hand, the combination Sx~ < —2, 
Sx't < —2 which corresponds to r\^{k~.^) < i, r\cr{k^.^) > 
^, cannot be realized for stable GSs. It is for these rea- 
sons that only for the two combinations of 6x~ and dx'^ 
explicitly presented in Eq. (153) can A~/AJ > 1 in prin- 
ciple amount to a non-trivial condition; below we shall 
demonstrate that indeed the latter cannot be trivially 
satisfied. For completeness, 5x~ = and Sx'^ = corre- 
spond to the general category of Sx~ > — 2 and (5x+ > — 2 
respectively for which, as indicated above, r\^{k~.^) > | 
and ncj{k'^.^) < i. In this connection, it is interesting 
to note that following Eqs. (146), (132), (122), (25), (26) 
and (17), Sx~ = 5x+ = corresponds to cases where 
(not exclusively) S"''(fe) is independent of fe (such as is 
the case for the uniform metallic GSs of the single-band 
Hubbard Hamiltonian), so that with reference to the find- 
ings in [6], the results r\a{k~.^) > \ and naik^.^) < i, 
obtained here, were in fact to be expected. The fact 
that, for 5xt = gcj/a<j — > one has ncr{k^.^) > \ and 
n<j(A;^^) < 5 implies that, for interaction potentials of 
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non-vanishing range, g/eo must be greater than some fi- 
nite value gc/eo in order for n<^(fc~^) < 5 or n<^(fc^^) > i 
to be realized (c/. footnote 15 above and footnote 23 be- 
low); here cq stands for an energy scale specific to the 
system under consideration. 

We now demonstrate that the condition A~/A+ > 1 in 
Eq. (153) cannot be trivially satisfied; this wc do through 
deducing restricted intervals in which various parameters 
of the problem at hand (such as a;^ and y^, and thereby 
6^, and d^) have to be confined in order for the men- 
tioned condition to be fulfilled. Readers who wish to take 
the validity of this statement for granted, may omit the 
remaining part of this Section. 

Making use of the property x~ > (see Eq. (135) and 
(II) in Eq. (127) above) we re- write the expression in 
Eq. (153) as follows 



where 



/(a;):=l + (1 + 4x(l + a; + ^a;))^/^ 



(154) 



(155) 



a-~ [1 - (1 + 4a;+(l + x+ + fe+)) ^'^\ . (156) 

X(j 

Above for clarity we have employed the short-hand no- 
tations f{x'^) and a for f{x~,Sx~) and a(a;+,i5a;+) re- 
spectively. When necessary, wc shall employ the latter 
more complete notations. We have 

f{x-)^2x- + {2 + 6x-)=:g{x-) for a;; ^ 00. (157) 

For the considerations that follow it is crucial to realize 
that in consequence of Eq. (137) and (II) in Eq. (127), 
x~ is bound to satisfy 



where 



:max{0, -(1 + (5.T^)}. (158) 



Consequently, unless indicated or implied otherwise, in 
what follows we implicitly assume that x~ satisfies the 
condition in Eq. (158). We have 



(159) 



The graphical representation of the functions on both 

sides of Eq. (154) clearly reveals the conditions under 
which the inequality in Eq. (154) (or Eq. (153)) is sat- 
isfied. One readily verifies that, for —2 < dx~ < 0, 
/(.T^) possesses a global minimum, equal to 1 -|- (l — 

(1 +_fe-)2)'/', at X- = -(1 + dx-)/2; for fe" < -2, 
f{x~) is complex- valued at x~ = —(1 + Sx~)/2. It is 
straightforwardly shown that for Sx~ < —2 and Sx~ > 0, 
f{x~) approaches the linear asymptote g{x~) (defined in 
Eq. (157) above) from below, whereas for —2 < 5x~ < 
it approaches g{x'^) from above. From these observations 
we deduce that the inequality in Eq. (154) is always vio- 
lated when the following conditions are satisfied: 



a>2 A Sx~ < -2. 



(160) 



On the other hand, the relationship in Eq. (154) is always 
satisfied when either 



a<2 ASx~ > -2 



or 



a < 



^{l + 6x^) 



A 5x~ < -2. 



(161) 



(162) 



Note that, for Sx' < -2, we have -2/(1 + Sx') e (0, 2]. 
One trivially obtains the following for a = a(a;+, (5a;+) 



'2 for Sx+< -2. 



(163) 



In this connection, note that in general (see Eq. (147) 
above) 



x+ < ~{l + 5xt\ 



(164) 



In the present case where a;+ < (owing to d+ < 
considered here; see Eq. (145) above), a;+ is bound to 
satisfy 



x+ <min{0,-(l-|-fe+)}. 



(165) 



Combining the above results, we conclude that the in- 
equality in Eq. (154) (or Eq. (153)) is violated when the 
following conditions are satisfied: 



fe;: < -2 A 5xt < -2. 



(166) 



These conditions are precisely those we excluded for sta- 
ble GSs in our discussions following Eq. (153) above. 

Considering the case corresponding to (c/. Eq. (153) 
above) 5.t^ > — 2 and bx'^ < —2 (the latter according 
to Eq. (163) implying a > 2), from our above consid- 
erations we deduce that the condition in Eq. (154) (or, 
equivalently, that in Eq. (153)) is satisfied provided that 
(see Eq. (158) above) 



x„ <x„ < 



2a + A{l + Sx-) 
a2-4 



a > 2. 



(167) 



As for the case corresponding to (c/. Eq. (153) above) 

5x~ < —2 and Sx'^ > —2 (the latter implying a < 2; 
see above, and in particular Eq. (163)), from our above 
considerations we similarly deduce that the condition in 
Eq. (154) (or, equivalently, that in Eq. (153)) is satisfied 
provided that 



-4{l + dx-) -2a 



-{l + 5Xa) 



<a<2. 



(168) 



Note that for a = -2/(1 fa") the RHS of the above 
inequality for expected, equal to — (1 -l- 6x^ ) = 

2 /a; the condition considered here, namely 6x~ < —2, 
implies that — (1 + Sx~) > 1. 
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6.2.2.2. The case when d+ > 

As we have indicated above, in the present case both 
solutions and A^^ of Eq. (129) are positive, with 
< Ajr_,_. For these solutions wo have 



A+±=£±(4,t/+), 



(169) 



where a;+ and y+ have been defined in Eq. (145) above. 
It is readily verified that 



A+ > 1 



A+ < 1 
■'V- > ^ 



Sx+ > - 2. 



(170) 



In view of the fact that Aj_ < AJ^_|_, it follows that we 
solely need to consider the case 

< A+_ < 1 < A++ for Sx+ > -2, (171) 

which in fact shows the necessity for 

c+ > (172) 

According to Eq. (142), for < -2o^/g, n„{k:^.^) < 
\ so that A^^ (which according to Eq. (152) is greater 
than unity for > —2ac/g and thus following Eq. (32) 
corresponds to nc{k'^.^) > ^) is not permitted. It follows 
that Ajr_ is the appropriate solution to Eq. (129) for A+; 
for this solution to correspond to a non- negative Z^^,^ , it 
is required that the following condition be satisfied: 



A 



^'^•^"'^"^ < 1 for 5x- < -2 A 5xt > -2. 



Aa C+{Xa,ya) 



(173) 



Here x^ naturally satisfies the inequality in Eq. (158) 
above. On the other hand, since in this Section we deal 
with d+ > 0, from Eqs. (145) and (147) we have (c/. 
Eq. (165) above) 



0<x+ < -{1 + Sx+). 



(174) 



The inequalities in Eqs. (171) and (174) restrict the range 
of variation in Sxt as follows: 



-2<fe+<-(l + a;+)<-l. 



(175) 



The condition A+./A" < 1 in Eq. (173) can be cast into 
the form f{x~) > ax~ presented in Eq. (154) above, 
with f{x) and a as defined in Eqs. (155) and (156) re- 
spectively. In spite of this, one should realize that here 
a;+ > (following d+ > 0; see Eq. (145) above), whereas 
in our previous considerations x+ < 0. 

As in the case of Eq. (153), below wc show that the 
condition A'^_/A~ < 1 in Eq. (173) cannot be trivially 
satisfied. Readers who wish to take the validity of this 



statement for granted, may omit the following, up to and 
including Eq. (179) below. 

Along the same lines as in the case of d+ < 0, we 
deduce that the inequality in Eq. (173) is always violated 
when 



a>2 A dx~ < -2. 



(176) 



On the other hand, the inequality in Eq. (173) is always 
satisfied when 



a < 



-{1 + SXa) 



A 6x~ < -2. 



(177) 



From the defining expression for a in Eq. (156) it can 
be readily deduced that for < a;+ < — (1 + (5a;+) (see 
Eq. (174) above), a can take non- negative values only for 

—2 < (5x-+ < —1. In fact we have 

< a < 2 for - 2 < 5a;+ < -1 A < x+ < -(1 + Sx+); 

(178) 

a = 2 applies only for Sx'^ — —2 and a; J € [0, 1]. Assum- 
ing that a;+ and fe+ are such that a > — 2/(1 -|- 6x~), 
it is readily deduced that the inequality in Eq. (173) is 
satisfied provided that 



x„ > 



-4(l + ^a;-) -2a 



-{l + 5xa) 



<a<2. 



(179) 



Following Eq. (142), for > —2aa/g (or, equivalently, 
5x~ > —2), ncr{k^.^) > i so that, unless Ajr^ > A~ 
(which entails Zk.^.^ < 0), Aj^_|_ cannot be a priori dis- 
carded. Evidently, since here A^_ > 0, in the event 
A < A~ , one will be confronted with a situation where 
there are two possible values to be taken by nrj{k^.^). 
Which of the two is the correct value for v\^{k^.^) is a 
matter that has to be decided on the basis of a further 
consistency test. Below we shall present an example of 
such a test but, before doing so, we derive the conditions 
to be satisfied for the occurrence of A < A^ . Readers 
who do not wish to consider this aspect, may omit the 
following, up to the last paragraph of this Section. 

It can be easily shown that for x~ > (which holds 
on account of (II) in Eq. (127) above) we have 

A+^<A- ^ f{x-)>px-, (180) 

where f{x) has been defined in Eq. (155) above and 



P:=—fixt,Sx+). 



(181) 



It is readily verified that, for a;+ > 0, /3 > 2. For a;+ > 
(which holds on account of our present assumption, 
namely (i+ > 0) one equivalently has 
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(T+ — O 



f{xt)<P'x+, 



where 



0':=—fix-,dx-). 

X(T 



(182) 



(183) 



It is easily deduced that for x~ > 0, 0' > 2. On account 
of 6x~ > —2 (see the text fohowing Eq. (179) above), 
along the lines of our earlier considerations we deduce 



f{x-)>Px-- 
Similarly 



^.^ 2^.4(l_Y».-) , ,,,,, 



2/3' + 4(1 + 6x+) 



/J'^ -4 



(185) 



Note that from the defining expressions in Eqs. (181) 

and (183) and the expressions in Eqs. (169) and (134) 
respectively (see also Eq. (133) above) wo have 



/? = 2A++, /3' = 2A;+ (^ 2A-). 



(186) 



In order to establish the conditions under which Aj_|_ < 
A~ is realized (if at all), it is natural first to establish the 
condition(s) required for the realization of Aj^ = A~. 
For A^_|_ = A~ to materialize, in which case (3 = (3' 
(following Eq. (186)), we must have 



2/3 + 4(1 + fa-) ^ ^+ _ 2(3 + 4(1 + 5xt) 



(3^-4 



(3^ 



■ (187) 



On the other hand, these expressions cannot be mean- 
ingful unless > — 1 (or x~ > x~] see Eqs. (137) and 
(158) above) and t/+ < —1 (sec Eq. (147) above), from 
which it follows that (note that (3 > 2) 

fe; > -1 - I A -2 < 5xt < -1 - |- (188) 

For /3 > 2, these inequalities arc more restrictive than 
5x~ > —2 and —2 < fa+ < —1 respectively. Note 
in particular that the range of variation in Sx'^ as re- 
quired for the satisfaction of Ajr^ = A^ becomes nar- 
rower, the closer to 2 that P approaches. The inequalities 
in Eq. (188) imply the following intervals for x~ and x+ 
corresponding to (3 — (3' (or A+_,_ = A^) 



2 2 
> - A 

- 13 13 + 2 



(189) 



One observes that in the event x~ = a;+, there is only 
one possibility for A^_|_ = A" to be realized, namely when 
x~ = x'^ = 2/(3. In view of the results in Eqs. (138) 
and (148), x~ = 3;+ in combination with A^^_ = A" 
(note that A^ = A;;^ since A;;_ < 0) imply A+_ = 
A~_ = 0. Since the choice of Ajr_ for A+ (which is the 
unrivalled choice in cases where d+ < 0) in the case when 



A^_ = implies that r\a{k'^.^) = 0, and this in turn 
amounts to the highly unusual condition for the single- 
particle spectral function, namely A„{k^.^;£) = for all 
e < fx, wc conclude that, in the event x~ = .xjjr and 
A+_,_ = A~, A+ is to be identified with A+_,_. If follows 
that, in such event, Zk^,^ = 0. Note that the inequalities 
in Eq. (189) show that x~ > .tJt which in consequence 
of the results in Eq. (187) implies that, for A'^_^_ = A~ to 
be realized, it is necessary that 6x~ > fa+; the equality 
x~ = .X implies the equality Sx~ = Sx'^ , which following 
the inequalities in Eq. (188) is only possible if 6x~ = 
Sx+ = 2/(3. 

It can be shown that f{x,Sx)/x (and therefore (3 and 
(3' for X = x^ ,6x = 6x'^ and x = x~ ,Sx = 6x~ , re- 
spectively; see Eqs. (181) and (183)) is a monotonically 
decreasing function of x (here for a; > 0) and a monoton- 
ically increasing function of 6x (here for 6x satisfying the 
appropriate inequality in Eq. (188) above). These prop- 
erties are properly reflected in the asymptotic expression 



—fix, 6x) ~ 2 + ^ for X 

x X 



00, 



(190) 



from which one further observes that so long as i5a; > —2, 
f{x, 5x)/x (and thus /3 and /3') never attains the value 2 
for any finite x. 

From the above properties it follows that, by start- 
ing from the set of values for x~ , Sx~ , x+ and Sx^ 
corresponding to Ajr_|_ = , one can achieve the in- 
equality Aj^_|_ < AgT by (if possible; see the inequalities 
in Eq. (189) above) increasing a;+ and/or (if possible) 
decreasing x(( while maintaining the values of and 
6x~ corresponding to Ajr_,_ = A~ (see Eq. (186) above 
according to which = 2A+^_ and /?' = 2A~). Alter- 
nativc^ly, AjJt^ < A~ can be achieved by maintaining the 
initial values of a;+ and x~ while (if possible) increasing 
5x~ and/or decreasing Sx'^ from their initial values. Ev- 
idently, these choices are not unique, and one can obtain 
(if at all possible) A < A~ by other apparent combina- 
tions of changes in the parameters corresponding to the 
condition A jr_^ — A~ . 

We now deal with a specific case which illustrates in 
how far the condition Aj_|_ < A^ can be considered as 
feasible; it turns out that for the specific case that we 
consider, Aj^_,_ < A~ is never realized. In § 6.3 below we 
investigate FL metallic states corresponding to two-body 
interaction potentials that are of shorter range than the 
Coulomb potential and arrive at two specific conditions 
(see Eqs. (195) and (196) below) which result in a;+ = 
— g6^/(aCTA+) and x~ = —gb^/{a„A~). From these we 
deduce that 

Sx- = -X- - A+X+, fe+ = -xt - A-X-, (191) 
which lead to 

5x- - 5xt = {A- - l)x- - (A+ - 1)4. (192) 
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From this it follows that, for A~ = A+ and in the event 
A+ were to be identified with A^T^ = j3/2, through em- 
ploying the expressions for x~ and a;+ in Eq. (187) one 
would have 



5x„ - 5xt = 



(3 + 2 



{6x^ - 6x+); 



this and the fact that /? > 2 imply that 

2 



5x„ = Sxt = -1 - 



/3' 



(193) 



(194) 



In arriving at the result on the RHS we have made use of 
the inequalities in Eq. (188) above, which uniquely deter- 
mine the only common value that Sx~ and 5x'^ can take 
in the case under consideration. Prom the expressions in 
Eq. (187) we observe that for 5x~ = Sx'^ = -1 - 2//?, 
indeed x~ = x+ = 2/(3 and consequently, through 
Eq. (191), A~ = A+ = (3/2. For completeness, here 
A^ = A;;_,_, A+ = A+^ and A~_ = A~_ = 0. We observe 
that, because in the present case x~ and x+ achieve re- 
spectively the highest and the lowest values that are com- 
patible with the inequalities in Eq. (189) (similarly and 
equivalently, for dx~ and according to the inequali- 
ties in Eq. (188) above), it is not possible to realize the 
strict inequality A+^ < A~ by appropriately altering the 
parameters (as described above) xf, Sx^ corresponding 
to A++ = A-. 



6.3. The case when 7 = 1, r = 1: Fermi Hquids 

As we have indicated earlier (§ 3), for metallic states 
of fermions interacting through the long-range Coulomb 
potential, "dk-.a is not diffcrcntiablc in the neighbourhood 
of Sp-a- For this reason, considerations in this Section 
should be viewed as not concerning the metallic GSs of 
the Coulomb-interacting fermion systems; the treatment 
appropriate to the latter systems is presented in § 6.1.3 
above. 

Along the lines of reasoning as in [6] concerning FL 
metallic states, we deduce that for such states we must 
have (compare with Eqs. (105) and (106) in [6]) 



bt 



-K4 



(195) 
(196) 



In the case of the uniform GSs of the Hubbard Hamilto- 
nian, where 6^ = ; the results in Eqs. (195) and (196) 
lead to A+A~ = 1 and consequently to n^ikT,-^) = \ 
for fepio- e S^-a (cf- Eq. (110) in [6]). For two-body in- 
teraction potentials which have longer range than the 
contact potential in the conventional Hubbard Hamil- 
tonian, Ca-{k^.^) cannot be identically vanishing for all 
kp-cr G Sp-a SO that, in general, &^ 7^ (sec Eqs. (123) 
- (125) above) whereby it is not necessary that for FL 
metallic states no-(fcF;<T) = ^ should hold. Further, in 



contrast with the case of the Hubbard on-site interaction 
for which under all conditions considered in [6] one has 
□^(fe^^) > 5 and n„{k'^.^) < ^, even for FL metallic 
states considered in this Section both ncr{k~.^) < | and 
r\„{k^.^) > i, satisfying, however, the necessary stabil- 
ity condition ncr{k^.^) — ncr{kp.^)=:Zk^.^ > 0, are feasi- 
ble. Below we explicitly demonstrate the feasibility of 

From Eqs. (125), (132) and (196) we obtain 



■d- 



(197) 



Prom (I) and (III) in Eq. (127) we have > -a^/g and 
6+ < —a^/g respectively so that we can write 

b- = -^,X-<l; 6+ = -^, A+>1. (198) 
g g 

On account of the expressions in Eqs. (197) and (198) it 
follows that, for 



Ar < 



2-A^ 



(199) 



one achieves the condition c~ < —2a^/g required for 
A~ < 1 (sec Eq. (142) above) and thus r\cr{k'^.^) < \. 
Evidently, the consistency of the condition A~ < 1 with 
that in Eq. (199) implies the further requirement that 



i< At < 2 - a: 



(200) 



which encompasses both A~ < 1 and A+ > 1 indicated 
in Eq. (198). It follows that, for two-body interaction 
potentials of non- vanishing range (corresponding to c~ ^ 
0), r\a{k'^.^) < ^ is feasible even for PL metallic GSs 
considered in this Section. 



§ 7. SUMMARY AND CONCLUDING REMARKS 

In this paper we have considered some exact properties 
of the uniform metallic GSs of systems of fermions for ar- 
bitrary spatial dimensions interacting through arbitrary 
isotropic pair potentials that are capable of being Fourier 
transformed; thus, for instance, systems of fermions in- 
teracting through the Lennard- Jones 6 — 12 potential lie 
outside the domain of applicability of the considerations 
in this paper, however those interacting through the Aziz 
et aJ. potential [104] lie inside. Our present work con- 
cerns natural extension of previous work [6] in which we 



Both of these potentials have been used for describing the 
inter-atomic interaction of ^He atoms, with the Aziz et al. 
potential as being the more superior; for a discussion of these 
and of a number of other available pair potentials see [104] 
and [105]. 
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have explicitly dealt with the conventional single-band 
Hubbard Hamiltonian where the interaction is operative 
only between particles of opposite spin at the same lattice 
site. 

In the present paper we have explicitly demonstrated 
that, for k approaching the Fermi surface <Sf;<t of the 
interacting metallic GS and s approaching the Fermi en- 
ergy Ep of the interacting system, the self-energy T,cr{k; e) 
continuously approaches the exact Hartree-Fock self- 
energy T,^^{k). On the basis of this and of an exact 
expression for the deviation of I]<j(fc;eF) from T,^^{k), 
obtained from the Kramers-Kronig relation expressing 
the real part of the self-energy in terms of the imaginary 
part, we have deduced the subsidiary equation Sa{k) = 
(Eqs. (52) and (51)), which in conjunction with the quasi- 
particle equation in terms of Yl^^{k) (Eq. (46)) fully de- 
termines Sp-cr (see Eq. (53)). Thus, whereas the latter- 
mentioned quasi-particle equation in terms of T,^^{k) is 
satisfied for all k on S-p-a, the reverse does not necessar- 
ily obtain; the failure to satisfy Sc,{k) = 0, in spite of 
satisfying the quasi-particle equation in terms of S™(fe), 
implies that for the k at hand the quasi-particle equation 
in terms of T,„{k;e) (Eq. (45)) is not fulfilled; follow- 
ing the arguments presented in [6], the set of all such k 
points constitutes the pseudogap regions of what other- 
wise would be constituent parts of the Fermi surface S^-^a- 
Amongst others, the GS momentum distribution function 
n„{k), though singular, is continuous for k transposed 
from inside to outside the underlying Fermi sea through 
pseudogap regions. For fe e St^-o- the expression corre- 
sponding to Sa{k) = amounts to a sum rule concern- 
ing Ini[Scr(fe;e)] (Eq. (54)); in this capacity, it can be of 
considerable relevance both for the purpose of both con- 
straining theoretical models for So-(A;; e) and determining 
the consistency of experimental results that are directly 
related with Tj„{k; e), or Im[So-(fe; e)] (such as those cor- 
responding to photoemission and invcrsc-photocmission 
experiments) . Owing to the generality of our considera- 
tions in this paper, the above-mentioned results equally 
apply to the uniform metallic GSs of the conventional 
single-band Hubbard Hamiltonian, a fact that can be ex- 
plicitly verified. 

In connection with the pseudogap phenomenon as ob- 
served in, for example, the angle-resolved photoemission 
data concerning the normal states of the cuprate com- 
pounds in the underdoped regime [106-110], it is relevant 
to mention that on general grounds it can be shown [78,6] 
(see in particular § 3.4 in [78]) that this phenomenon can- 
not be associated with the uniform GSs of systems for 
which k is defined over the unbounded reciprocal space; 

thus the single-band Hubbard Hamiltonian (for which 
the reciprocal space is confined to the IBZ of the under- 



In these cases, indirect gap is not excluded. Further, for 



lying Bravais lattice), even in cases where the interaction 
is operative solely amongst particles of opposite spin on 
the same lattice site, is special in that, through the pe- 
riodicity over the IBZ of the corresponding functions of 
k, it takes account of Umklapp processes, albeit in some 
restricted form; this aspect is vital for the possibility of a 
uniform metallic GS of the single-band Hubbard Hamil- 
tonian supporting the pseudogap phenomenon. With ref- 
erence to the observations in [6], for such states <Sf|2\>5f;(7 
is non-empty and, in view of our observations in the 
present work, S„{k) ^ for all k G »SF°i\<SF;(7- 

For the uniform and isotropic GSs of the uniform- 
electron-gas system, where on account of isotropy and of 
the Luttinger theorem concerning the 'volume' enclosed 
by <Sf;<7 [92,80] one has Sj^-a- = "Sp-l^ through making use 
of S(^(fe; s) ^ S^''(fe) for fc approaching S^^^ and s £f, 
we have deduced (§5) exact expressions for the exchange 
potential /Ltx;cr and correlation potential /Uc;<t) contributing 
to the exchange-correlation potential Atxc;(7 = Mx;(t + A*c;(t 
which plays a central role within the framework of the 
GS density- functional theory [75-77]. Our expression 
for /Ux;<T) specialized to systems of electrons interacting 
through the Coulomb potential for = 3, identically re- 
produces the well-known exact result for the exchange 
potential. Making use of the RPA for n<^(fe) correspond- 
ing to the paramagnetic GSs of the uniform-electron-gas 
system (note that, for systems of Coulomb-interacting 
fermions, n^'''^(fe) — n'^\k) reproduces the leading-order 
contribution to the exact no-(fe) — nlj\k) in the high- 
density limit), we have numerically calculated Hc;cr which 
in the high-density limit, and up to densities correspond- 
ing to the dimensionless Wigner-Seitz radius r^^l, is 
in excellent agreement with the /ic;cr deduced from the 
quantum Monte Carlo results for the GS correlation en- 
ergy (see Figs. 1 and 2). 

We have explicitly shown that interactions of non-zero 
range are of non-trivial consequence to a number of GS 
and excited-states properties of metallic systems. For 
instance, whereas for a strictly contact-type two-body 
interaction potential ncr{kp.^) > j and n^{k'l.^) < i 
for feF;(7 € <SF;cr [6], for two-body interaction potentials 
of non- vanishing range in principle both n^rik^.^) < i 
and n„{k^.^) > i are feasible (in both cases, subject to 
the condition r\r^{kp.^) — r\a{k'^.„) > 0, implied by the 
assumed stability of the underlying GSs). This aspect 
is borne out by the quantum Monte Carlo results for 



ferromagnetically ordered uniform GSs one encounters direct 
spin gap even in models where k is defined over the unbounded 
reciprocal space. It is interesting to note that the approximate 
solution of the Hubbard model, the so-called Hubbard-I solu- 
tion [10], predicts an indirect (Mott-Hubbard) gap (see Fig. 1 
in [10]). 
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no-(fe) pertaining to the uniform and isotropic GSs of the 
Coulomb-interacting electron-gas system for d = 2 [111] 
(see Fig. 5.9 herein) in the low-density regime and of 
^Hc in the normal liquid state at sufficiently large am- 
bient pressures [105,113-116]. In the light of this and 
of the strict inequality ncr{k'^-a-) > 5 applicable to the 
na-{k) pertaining to the uniform GSs of the conventional 
single-band Hubbard Hamiltonian [6] , employed in stud- 
ies of ^He [117-119], as well as to the uniform GSs of the 
Stoncr Hamiltonian (involving contact interaction but 
defined on the continuum; likewise employed in studies 
of ^He within the framework of the paramagnon theory 
[120-124]; for a review of spin-fluctuation theories sec 
[125]), one observes that such Hamiltonians cannot be 
capable of reliably describing the normal liquid state of 
^He. In this connection, it is in place to refer to a series of 
papers by Hirsch (see [55] and references herein) in which 
the viewpoint has been consistently advocated that many 
(collective) phenomena for whose rationalization tradi- 
tionally the conventional single-band Hubbard Hamilto- 
nian has been invoked should in fact be accounted for 
by means of 'extended' Hubbard Hamiltonians (such as 
considered in [49-65], or in [126-129]), involving site 
off-diagonal (to be contrasted with on-site) interaction 
terms (including, e.g., the nearest-neighbour exchange). 
In [55], Hirsch referred to such phenomena as itinerant 
ferromagnetism [130,131] (see also [6]), heavy-fermion 
superconductivity [132,133], high-temperature supercon- 



We recall tfiat, for d = 3 and Coulomb- interacting 
fermions, n'^''^(kp.,y) submerges below i for densities corre- 
sponding to Ts > 6.09887... (see footnote 15 above); the 
Monte Carlo results by Ortiz and Ballone [112] on the other 
hand show that, for rs = 10, ncr{kp.a-) ~ 0.69. Interestingly, 
for d = 2, according to the quantum Monte Carlo results 
by Conti [111], for rs — 10 one has nc^{k^.,y) ~ 0.55 and for 
Vs = 20, na{k^;a) ~ 0.42. We observe that, for sufficiently 
large interaction strength, indeed the non-vanishing range of 
the interaction potential (here the Coulomb potential) gives 
rise to no-(fe^;,-) < 5- 

Concerning the role of interactions of longer range than 
the on-site repulsion in promoting specifically partial ferro- 
magnetism, we refer the reader to the phase diagram in Fig. 4 
of [62] which illustrates, albeit on the basis of a mean-field ap- 
proximation, the consequence of a nearest-neighbour repuJsive 
interaction energy V (note that in the above-mentioned phase 
diagram the U axis intersects the V axis at some positive 
value of V) . In this connection we mention that according to 
the analysis in [6] , in the case of the conventional single-band 
Hubbard Hamiltonian, partially polarized uniform ferromag- 
netic states are often barred from qualifying as eigenstates 
through a kinematic constraint (sec Eq. (63) in [6]); it can 
be readily verified that for longer-range interaction potentials 
no such a priori kinematic constraint exists, this in conse- 
quence of the k dependence of S^''(fc) in the case of the latter 
potentials. 



ductivity in the cupratc compounds [44] (see also [134]), 
and superfluidity in hquid "^He [55]. 

Since for the (correlated) GS kinetic energy we have 
^ eerier (fc), by continuity it follows that in general 
the GS kinetic energy of systems in which r\cr{k^.„) < \, 
kp-a £ Sf-ct, exceeds that of comparable systems in which 
rtaik^.a) > h (^ot*^ that n^(k) > and n^ik) = N, 

independent of interaction). We thus observe a direct 
relationship between the range of the two-body interac- 
tion potential and the magnitude of the corresponding 
GS kinetic energy, with the latter tending to increase (in 
the manner specified above) in consequence of increasing 
the range of the interaction potential. The existence of 
such a direct relationship between the range of the inter- 
action potential and the value of the GS kinetic energy 
should be of relevance to a better understanding of the 
physics of the cuprate superconductors for which there 



We do not unreservedly subscribe to the assignments 
in [55] quoted above. For instance, in our opinion it is 
not the phenomenon of superconductivity in heavy-fermion 
compounds that has been rationafized by means of the con- 
ventional Hubbard model (explicitly, in [135], to which in 
[55] is referred, one does not explicitly rely on the Hubbard 
Hamiltonian) . Rather, in considering the heavy-fermion com- 
pounds and superconductivity [132,133,136,137,135] in these, 
one has in essence solely relied on the so-called 'almost- 
localized Fermi-liquid' picture [117,118] of the charge carriers; 
the latter has its root in the Brinkman-Rice [138] scenario 
coucc^rning the nature of the metallic states of the Hubbard 
Hamiltonian (for d = 3) at half-filling for the on-site U ap- 
proaching from below a critical value Uc, for which, accord- 
ing to the Gutzwiller Ansatz for the GS wavefunction of the 
Hubbard Hamiltonian combined with the Gutzwiller approx- 
imation [9,139], the jump in the GS momentum distribution 
function at iSf;^ vanishes and thereby the inverse of the quasi- 
particle effective mass (hence the qualification 'almost local- 
ized'). For completeness, we mention that the Gutzwiller ap- 
proximation (to be distinguished from the Gutzwiller Ansatz) 
is a well-known source of uncertainty. For instance, the ex- 
act results corresponding to the Gutzwiller Ansatz for the 
GS wavefunction of the single-band Hubbard Hamiltonian for 
d = 1 [140] reveal that at half-filling the above-mentioned 
jump in the momentum distribution function does not van- 
ish at any finite value of U; in contrast, according to the 
Gutzwiller approximation for d = 1, this jump vanishes for 
Uc/t « 10 (see Fig. 13 in [140]); these results should in turn 
be contrasted with the exact result indicated in footnote 20 
above. Finally, insofar as the normal states of the heavy- 
fermion compounds are concerned, although the considera- 
tions in [141] are based on the Gutzwiller Ansatz (together 
with the Gutzwiller approximation) for the GS wavefunc- 
tion of the Hubbard Hamiltonian, the work in [142], based 
on the periodic Anderson Hamiltonian, leaves no doubt that 
the analogies proposed in [141] have indeed been based on 
"intuitive grounds" . 
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is growing evidence indicating that in the normal state 
the kinetic-energy contribution to the GS total energy is 
relatively larger than in conventional metals and that su- 
perconductivity in the former is driven by a consequent 
decrease in the kinetic energy of these systems upon en- 
tering into the superconducting state. In the light of this 
and of our theoretical finding, wc believe that any theo- 
retical modelling of in particular the cuprate compounds 
is required to take account of the non-zero range of the 
interparticle interaction potential. 

In addition to establishing the above-mentioned im- 
pact that a non-zero range of interaction has on the be- 
haviour of n^{k) in the vicinity of Sp-a, we have explic- 
itly exposed (§ 6) some distinctive features associated 
with the actual range of the interaction potential; that 
is whether this range is sharply limited to a finite inter- 
particle distance or it extents to arbitrary large distances. 
For instance, for interaction potentials v{r — r') that pos- 
sess power-law decay for \\r — r'|| oo, the singular be- 
haviour of ncr(fc) in the neighbourhood of any fe directly 
implies a singular behaviour in the Fock part ^^{k) of 
the self-energy in the neighbourhood of the same fe; the 
two singularities are, however, of different types. This 
observation is specifically significant for the fe points in 
the neighbourhood of iSf;<t, where ria^k) is uncondition- 
ally singular (and not necessarily discontinuous). From 
our considerations it follows that, for systems of fermions 
interacting through a potential v(r — r') decaying, for 
instance, like l/||r — r'|l"' for \\r — r'\\ ^ oo, a discon- 
tinuity in r\a{k) at a general fe = feF;<7 € >Sf;<7 implies a 
logarithmic divergence at k^-a- in the mth derivative (in 
the direction of the radius vector feF;^) of S^(fe). 

We believe that a detailed examination of r\cr{k) deter- 
mined experimentally (for fe in the close vicinity of iSf:<j), 
in conjunction with the appropriate expressions for this 
function as presented in this paper, can shed light on 
some essential aspects of the two-body interaction poten- 
tial to be employed within the framework of an extended 
single-band Hubbard Hamiltonian. The details presented 
in this paper can also with advantage be employed to con- 
strain (approximate) theoretical frameworks concerning 
correlation functions (such as the self-energy Scr(fe;e)) 
pertaining to the metallic GSs of the many-body Hamil- 
tonians considered in this paper. 
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APPENDIX A: ESTIMATION OF ERROR IN 
S^{k) WITHIN THE FRAMEWORK OF A 
FINITE-ORDER MANY-BODY PERTURBATION 
THEORY FOR E^(fc;e) 

In the defining expression for Sa{k) in Eq. (51) we 
encounter lm.[T,a-{k;SF ± s)] in which for definiteness we 
assume T,„{k;e), e e (—00,00), to have been calculated 
within the framework of a finite-order perturbation the- 
ory based on skeleton self-energy diagrams (for the def- 
inition see [92]) evaluated in terms of the exact single- 
particle Green fimction. With the latter assumption we 
wish to estimate the error in S„{k) as arising from the 
truncation of the set of skeleton self-energy diagrams and 
not from the further approximation corresponding to us- 
ing an approximate single-particle Green function. 

We proceed by writing S„{k) as follows 

5,(fe) = SW(fe;A) + S(2)(fe;A), (Al) 

where 

5W(fe;A):=_ / — Im[S<,(fe;eF + e) + I]^(fe;£F-e)], 

Jo £ 

(A2) 

1 dp 

5^2)(fe;A):=- / — Im[I]^(fe;£F+£) + Sa(fe;£F-e)], 

(A3) 

in which A > stands for an energy sufficiently large 
with respect to a relevant energy scale in the system 
under consideration [78] so that a finite-order large- |e| 
asymptotic series for Im[Scr(fe;eF i e)] is accurate for 
\e\ > A. Later in our considerations we shall consider 
the case in which the latter series is formally summed to 
all orders whereby it is possible to effect the limit A | 0. 
To avoid unnecessary technical complication, in what fol- 
lows we assume that the two-body interaction potential 
f(r — r') is of shorter range than the Coulomb potential 
and moreover that the single-particle excitation spectrum 
£fe is unbounded from above so that the single-particle 
spectral function Acr{k;s) is of unbounded support. For 
such cases we have [78] 

Im[S,(fe;£)]~=^ + ^^|^ + ... for \e\ ^ 00, 

(A4) 

where 



For the stability of the GS of the system under consider- 
ation it is required that S^{k) > 0, Vfc,a [78]. 



24 



E„{ky.=— / de' {Re[S<,(fe; s') + S<,(fe; -s')] 

-2S,;ooo(fc)}, (A5) 

-1 f°° 

U^{k):=— / de' {e'Rc[E„(fc; e') - -e')] 



(A6) 



In Eqs. (A5) and (A6) the real-valued functions 
S(j.oo^(A;), m = 0,1, are the coefficients of the ffist 
two leading terms in the large- 1 2; | asymptotic series of 
SctC^; z), with Im(2;) ^ 0, that is 



\z\ — > 00. 



(A7) 



Here Scr(fc;z) stands for the analytic continuation of 
Scr(fe; e) into the physical Riemann sheet of the complex 
z plane, for which we have T,„{k; s) = lim^^o ^a{f^> si**?)) 
e^H- We point out that, even for cases where T,„.^^{k) 
is bounded for all finite values of m, use of a complex z in 
Eq. (A7) is necessary when the series on the RHS of this 
equation is taken account of to some infinite order, be- 
cause Eo-(fc; z) possesses a branch-cut discontinuity along 
the real e axis. It can be shown that [82,78] 



We write 



S<,(fe)=sW(fe;A) + 42)(fe.A), 

n<,(fe) = n«(fc;A) + ni2)(fe;A), 



(A8) 



(A9) 
(AlO) 



where E^^\k;A), E^^\k;A), ui^\k;A) and Il^a\k;A) 
stand in relation to ^^(fc) and Il^{k), defined in 
Eqs. (A5) and (A6), as Si^\k) and si^\k) in Eqs. (A2) 
and (A3) in relation to Sa{k), defined in Eq. (51). Prom 
the defining expressions for the functions on the RHSs 
of Eqs. (A9) and (AlO), making use of the asymptotic 
series for Ti^ik; z) in Eq. (A7), for A ^ 00 one straight- 
forwardly obtains 



ttA 



37rA3 



TT(2)('fc. A^ ^^a;oo3(^) _ 2So-;oo5(fc) _ 

Prom the above results we deduce 

Im[Eo-(fe; £f + e) + Scr(fe; Ep - e)] 
nW(fc;A)-SpSW(fe;A) 

+ (S(j;o03 (fe) ~ £f 5]ct;002 C^)) 

\e\, A — > 00. 



(All) 
(A12) 



(A13) 



Making use of this expression, from the defining expres- 
sion for S''a '{k) in Eq. (A3) we obtain 

^(2)(^. ^) ^ [nW(fc;A)-£pSW(fc;A)]/7r 

+ ^[^— (^)]/-%...,A^oo. (A14) 

It is interesting to note the occurrence in Eq. (A14) (in 
particular in the second term on the RHS) of £f, the en- 
ergy with respect to which the energies of the lowest-lying 
single-particle excitations arc measured, in conjunction 
with coefficients that determine the behaviour of Ti„(k\ e) 
at large values of |c|. 

The expression on the RHS of Eq. (A14) is not a well- 
ordered asymptotic series, as the first term involves con- 
tributions to terms of higher order in 1/A than 1/A^. 
Since {li^}\k- A) - n,(fc)) = o(l) and (si^^(A;:A) - 
So-(fc)) = 0(1) for A — > 00, to the leading-order in 1/A 
we have 



5(2) (fc; A) 



[n,(fc) -£pS<,(fc)]/7r 



A2 



for A — »■ 00. 



(A15) 



Note that by replacing Wa{k) and ^^(A;) on the RHS 
of this expression by the expressions on the RHSs of 
Eqs. (A9) and (AlO) respectively and using the leading- 
order terms in the asymptotic series on the RHSs of 
Eqs. (All) and (A12) respectively, we recover the ex- 
pression on the RHS of Eq. (A14). 

Prom our perspective, the most important aspcx't of the 
asymptotic expression in Eq. (A14) is the explicit depen- 
dence on I]o.;oo2(fc) and T^^-ryo^ik) of the second term on 
the RHS. This in view of the fact that according to the 
analysis in [78] , the reproduction of the exact So-;oo2 (fe) 
requires calculation of Eo-(fe;e) to third order and that 
of Ti^-ao^ik) to fourth order in the perturbation series 
in terms of skeleton self-energy diagrams, with these in 
turn evaluated in terms of the bare particle-particle inter- 
action function and the exact interacting single-particle 
Green function. Thus the second term on the RHS of 
Eq. (A14) associated with the Scr(A;; e) evaluated in terms 
of the full set of first- and second-order skeleton dia- 
grams deviates from its exact counterpart. In connec- 
tion with the result in Eq. (A15) (or Eq. (A14)), it is 
important to realize that since s'^\k] A) and s!^\k; A) 
as defined in Eqs. (A2) and (A3) respectively exist for 
A J, 0, and since s!j'\k; A) ^ for A J, 0, one can 
in principle recover the entire S„{k) by evaluating the 
largc-A asymptotic scries of S^\k; A) to all orders in 
1/A and, subsequent to its full summation, taking the 
limit A J, 0. Although this approach is of no practi- 
cal use, the fact that the formal infinite-order large-A 
asymptotic series of ^^^^(fe; A) is capable of exactly re- 
producing Sa{k) is significant in that it sheds light on 
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the way (as in Eq. (A14) above) in which the coefficients 
{SCT;oo„(fe)} of the large-|£| asymptotic series of Y,cr{k;e) 
contribute to those of the former asymptotic series; by 
the further knowledge that for m > 1 all skeleton self- 
energy diagrams from order 2 up to and including order 
m + 1 contribute to S(^;oo„(fe) [78], one gains insight into 
the consequences for Sa{k) of employing a finite-order 
perturbation series for ^^(fc;^). The same line of rea- 
soning applies to (fe; A) and 11^ {k; A) from whose 
complete large- A asymptotic series the exact ^^(fe) and 
Ila{k) respectively can be deduced; these asymptotic se- 
ries are evidently fully determined by {T,a-oomi^)} (see 
Eqs. (All) and (A12) above). In this connection note 
that according to Eq. (A15) the leading-order term in the 
large-A asymptotic series of Sa {k) is in addition to 
fully determined by E^{k) and II„{k). Concerning e^, in 
view of Eqs. (46) and (A8) one observes that this constant 
is determined by the first-order (skeleton) self-energy di- 
agram; its exact determination is, however, dependent 
on the knowledge of the exact interacting single-particle 
Green function Ga(k;£) or, what in this context is the 
same, no-(fe). 

The inaccuracy in a calculated S„{k), which in prac- 
tice arises from use of both a finite set of (skeleton) 
self-energy diagrams and an approximate single-particle 
Green function (either a non-interacting single-particle 
Green function or one determined self-consistently), com- 
bined with the inaccuracy in S^^(A;) arising from use of 
an approximate single-particle Green function (or nc{k); 
see Eqs. (12), (16) and (17)), leads through Eqs. (45) 
and (50) to deviation of a calculated Fermi surface from 
its exact counterpart Sp-a- For the uniform paramag- 
netic metallic GSs of the single-band Hubbard Hamilto- 
nian where T,™{k) is equal to the exactly-known constant 
value Un/{2h), in which n = N/N^^ (the total number of 
fermions per lattice site), the inaccuracy in a calculated 
Fermi siufacc based on the use of Eq. (45) is entirely at- 
tributable to that in the underlying Sa{k) (see Eq. (50)). 
The apparent violation (albeit very slightly) of the exact 
result iSf;o- C S^l, deduced in [6], by the S-p-a (pertain- 
ing to the uniform paramagnetic GS of the single-particle 
Hubbard Hamiltonian, with U /t — 4) determined from a 
self-consistently evaluated T,„{k] e) (based on the second- 
order expansion of E<7(fe; e) [143]) is thus attributable to 
the fact that none of Scr;oo„ , m>2, corresponding to a 
second-order perturbation expansion of So-(fc;e) can be 
correct [78] . For our criticism of non-selfconsistently cal- 
culated Fermi surfaces pertaining to the metallic states of 
the single-band Hubbard Hamiltonian we refer the reader 
to [6] . Note the principal advantage of the expression for 
<Sf;o- in Eq. (53) in comparison with that in Eq. (45); 
following the evident fact that Y,^{k) is in principle ex- 
actly determined by the first-order self-energy operator 
(i.e. leaving aside the necessity for the knowledge of the 
exact no-(fe)), use of the expression in Eq. (53) in con- 



junction with a first-order expression for the self-energy 
operator guarantees the exact reproduction of the skele- 
ton of the underlying Fermi surface (if not the entire 
Fermi surface). □ 
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FIG. 1. The correlation potential fl^-cr (in Hartree 
atomic units) pertaining to the paramagnetic GS of the 
Coulomb-interacting uniform-electron-gas system for d = 3 as 
calculated through the expression in Eq. (87) with ncr{kFx) 
herein replaced by n'^''^{kpx), due to Daniel and Vosko [95], 
and the flcia determined through /ic;o- = — (rs/3)d£c/drs 
(c/. Eq. (92)) with herein replaced by a parametrized ex- 
pression due to Vosko ct aJ. [97] based on the quantum Monte 
Carlo results (curve QMC) due to Ceperley and Adier [99]. The 
open circles indicate the explicitly calculated data and the lines 
through these are cubic-splines interpolations. Explicit calcu- 
lation shows that ~ 0.0310907 In(rs) - 0.0570205 for 
rs J, 0. A least-squares fit of the expression aln(rs) -I- fe to 
curve QMC over the range presented yields a — 0.0308213, 
b = —0.0589913; a similar least-squares fit to curve RPA over 
the range presented yields a = 0.0304144, b = -0.170157. Ex- 
plicit analysis of the eight presented data points for p,^^^ and 
yields p?™ - /i^™ = 0.107 ± 0.002, to be compared with 
6|qmc — &|rpa ~ 0.111. The constancy of /Xc-CT^ — Mcf(T over the en- 
tire range presented is consistent with the fact that the potential 
/ic;(T is determined up to an additive constant. For completeness, 
a somewhat more accurate, but vastly time-consuming, calcula- 
tion of /i[!f* at four different rs values inside [10"'', 10"''] reveals 
that while p.^^^ — fll':^ remains independent of rs, the absolute 
value of 6|rpa decreases slightly (by about 0.003) with respect to 
the value presented above. For this reason, we do not consider 
6|rpa = —0.170157 as being the converged value. 
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FIG. 2. The correlation potentials /lefCT and ^5?;"*^ over a 
larger range than in Fig. 1 (for details of the underlying cal- 
culations see the caption of Fig. 1). The vertical broken line 
at ln(rs = 7.769269 . . .) = 2.050176 . . . indicates the boundary 
beyond which = x^J^{k'^.„) — 'n^^'^ (kt,a) is negative, in- 

dicative of the instability of the uniform paramagnetic GS of the 
system under consideration according to the RPA (see footnote 
15). Note the remarkable near-constant behaviour of /i[?™ — 
over the range of validity of the RPA. That this aspect cannot 
be due to smallness of higher-order correlation effects is evi- 
dent from the behaviour of finite-order asymptotic series of /ic;cr, 
for ra \ 0, at large values of rs, for which we employ (bro- 
ken curves) p^.^ ~ 0.0310907 In(rs) - 0.0570205 for curve (a), 
/ic;<T ~ 0.0310907 ln(rs)-0.0570205-0.0585 rs for curve (b) and 
Me;,, ~ 0.0310907 ln(rs)-0.0570205-0.0585 rs +0.2794 rs^^ for 
curve (c). These asymptotic series correspond to /ic™ . Qur cal- 
culation of /ic;o- based on the expression in Eq. (87) with naiknx) 
therein replaced by a parametrized expression due to Gori-Giorgi 
and Ziesche [100] yields values between approximately —5 x 10"^ 
and 7 x 10"'^ over the range rs ^ 12 where this expression is as- 
serted reliably to represent the momentum distribution function. 
The reason underlying this shortcoming is discussed in the text. 
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FIG. 3. The behaviours of a^:=[l — ncrik^^a-)]/^ (curve (a)) 
and ■.=ncr{k^^„)/a (curve (b)), where a~rs/(7r^7o), in which 
7o is defined in Eq. (88); the results presented here correspond 
to na{k) = n^™(fc) [95]. One observes that, although for suffi- 
ciently small Vs, Oct and approach towards the same value Qo-, 
in contrast with the finding by Daniel and Vosko [95] this value 
is not equal to approximately 1.7, a value to which and 
are indeed relatively close for 10~*^ Tb ^ 1. From the results de- 
picted in this figure, one further observes that, for Ts J, (explic- 
itly, for Ts H 10"'') n^''^ (ky.cr) , as defined according to Eq. (106), 
approaches | more rapidly than would be expected from the 
above-mentioned finding in [95]. 
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